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Abstract. The periodic Anderson model (PAM) captures the essential physics of heavy fermion materials. 
Yet even for the paramagnetic metallic phase, a practicable many-body theory that can simultaneously 
handle all energy scales while respecting the dictates of Fermi liquid theory at low energies, and all in- 
teraction strengths from the strongly correlated Kondo lattice through to weak coupling, has remained 
quite elusive. Aspects of this problem are considered in the present paper where a non-perturbative local 
moment approach (LMA) to single-particle dynamics of the asymmetric PAM is developed within the 
general framework of dynamical mean-field theory. All interaction strengths and energy scales are encom- 
passed, although our natural focus is the Kondo lattice regime of essentially localized /-spins but general 
conduction band filling, characterised by an exponentially small lattice coherence scale wl- Particular em- 
phasis is given to the resultant universal scaling behaviour of dynamics in the Kondo lattice regime as an 
entire function of u/ = lj/ujl , including its dependence on conduction band filling, /-level asymmetry and 
lattice type. A rich description arises, encompassing both coherent Fermi liquid behaviour at low-a/ and 
the crossover to effective single-impurity scaling physics at higher energies — but still in the o;/o;_L-scaling 
regime, and as such incompatible with the presence of two-scale 'exhaustion' physics, which is likewise 
discussed. 



PACS. 71.27.+a Strongly correlated electron systems; heavy fermions 
pounds and alloys; Kondo effect, valence fluctuations, heavy fermions 



75.20.Hr Local moment in com- 



1 Introduction. 

The paramagnetic metallic phase of heavy fermion ma- 
terials provides a classic example of strongly correlated 
electron physics [1,2]. Spin-flip scattering of itinerant con- 
duction electrons by essentially localized /-level electrons 
leads to large effective masses and the low-energy scale(s) 
symptomatic of any strongly correlated state. At low ener- 
gies and/or temperatures the lattice coherence is paramount 
and the system is a Fermi liquid with well defined quasi- 
particles and coherent screening of the /-level spins; be- 
haviour that crosses over for sufficiently high energies to 
essentially incoherent screening and the effective single- 
impurity characteristics of the Kondo effect [1,2]. 

Handling theoretically the many sides and attendant 
issues of this basic physics is of course another matter. The 
paradigm here is the periodic Anderson model (PAM), the 
natural lattice generalization of the Anderson impurity 
model (AIM), in which each lattice site contains a corre- 
lated /-level (with interaction U) hybridizing locally with 
a non- interacting conduction band [1,2]; and a description 
of which remains a major challenge, particularly in the 
strongly correlated Kondo lattice regime of effectively lo- 
calized /-spins but arbitrary conduction band filling. That 
reflects in large part the inherent difficulties in developing 
a many-body theory that can capture non-perturbatively 
the strong coupling regime of primary interest, satisfying 



in particular the dictates of Fermi liquid theory at low 
energies and yet capable of describing the problem on all 
energy scales. Moreover, no matter how strong the correla- 
tions, the Fermi liquid nature of the ground state implies 
adiabatic continuity to the non-interacting limit; so the 
same theory should also be able to handle the full range 
of interaction strengths, including simple perturbative be- 
haviour in weak coupling. 

Our aim in the present paper is to develop an ap- 
proach to the paramagnetic phase of the PAM that meets 
the above criteria, within the general framework of dy- 
namical mean-field theory (DMFT) [3-6]. The PAM has 
of course been studied extensively within DMFT using 
an impressive range of methods. Numerical techniques in- 
clude the numerical renormalization group (NRG) [7,8] 
quantum Monte Carlo (QMC) [9-12] and exact diagonal- 
ization [13]. Theoretical approaches range from perturba- 
tion theory in the interaction U [14,15], iterated pertur- 
bation theory [16,17], the lattice non-crossing approxima- 
tion [18,19] and the average i-matrix approximation [20], 
to large- N mean-field theory/slave bosons [21-23] and the 
Gutzwiller variational approach [24,25]. Such techniques 
nonetheless possess well known limitations [2]; be it an 
inability to handle strong correlations, failure to recover 
Fermi liquid behaviour or even the non-interacting limit, 
unrealistic confinement to the lowest energy scales and so 
on. NRG aside, analogous comments apply to full scale nu- 
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mcrical methods. QMC for example is restricted to mod- 
est interactions and relatively high temperatures, while 
finite-size effects render exact diagonalization of limited 
value. These remarks are certainly not intended to de- 
tract from the many insights that have accrued from such 
approaches. They are made simply to emphasise the desir- 
ability of developing new, necessarily approximate theories 
for this longstanding problem. 

One such is pursued here, the local moment approach 
(LMA) [26-35] , the primary emphasis of which is on single- 
particle dynamics and transport. Initially developed in the 
context of pure quantum impurity models (AIMs) [26-33], 
the LMA is intrinsically non-perturbative but technically 
quite simple, with the physically intuitive notion of lo- 
cal moments introduced explicitly from the outset. This 
leads naturally to an underlying 'two-self-energy' descrip- 
tion in which the essential correlated spin-flip physics is 
readily captured; and corresponds physically to dynamical 
tunneling between initially degenerate local moment con- 
figurations, which in lifting the spin-degeneracy restores 
the local singlet symmetry characteristic of the Fermi liq- 
uid state. The desiderata mentioned above are well met 
[26-33], all interaction strengths and energy scales be- 
ing encompassed, including the low-energy requirements 
of Fermi liquid theory (although the approach can also 
handle models with non-Fermi liquid phases, see e.g. [31- 
33]). In particular, for the strong coupling Kondo regime of 
the conventional metallic AIM, LMA results for dynamics 
have been shown [28,29,33] to give very good agreement 
with NRG calculations; and, for static magnetic proper- 
ties, with exact results from the Bcthe ansatz [30]. 

More recently, exploiting the fact that within DMFT 
all correlated lattice-fermion models reduce to an effec- 
tive quantum impurity hybridizing self-consistently with 
the surrounding fermionic bath [3-6], the LMA has been 
extended to encompass the particle-hole symmetric PAM 
[34,35]. Here the system is ubiquitously a 'Fermi liquid in- 
sulator' that evolves continuously with increasing interac- 
tion strength from a simple non-interacting hybridization- 
gap insulator to the strongly correlated Kondo insulating 
state; with an insulating gap scale that becomes exponen- 
tially small in strong coupling, such that physical prop- 
erties exhibit universal scaling in terms of it (i.e. contain 
no explicit dependence on the 'bare' high-energy material 
parameters, U etc, that enter the underlying model Hamil- 
tonian) . A comprehensive description of single-particle dy- 
namics [34,35], electrical transport and optical properties 
[35] of Kondo insulators arises, encompassing all relevant 
frequency (u>) and/or temperature (T) domains; and ex- 
ploitation of scaling in particular enables direct, rather 
successful comparison to a range of experiments [35] . 

Important though it is to the problem of Kondo insu- 
lators the particle-hole symmetric PAM is of course spe- 
cial, and the desirability of developing the LMA to en- 
compass the asymmetric PAM and hence the generic case 
of heavy fermion metals is self-evident. That is consid- 
ered here, our specific focus being on T = single-particle 
dynamics. In addition to intrinsic interest in such per se, 
and the fact that their u-dependence exemplifies much 



of the underlying physics of the problem, knowledge of 
single-particle dynamics is well known [3-6] to be suffi- 
cient within DMFT to determine q = transport and op- 
tical properties; which will be considered in a subsequent 
paper (in parallel to previous LMA work on Kondo in- 
sulators [34,35]). The present paper is accordingly organ- 
ised as follows. After appropriate background to the PAM 
within DMFT (§2), formulated for an essentially arbitrary 
lattice, implications of adiabatic continuity and the Lut- 
tinger integral theorem [36] are considered in §3; together 
with the quasiparticle forms for the local conduction (o) 
and /-electron spectra that Fermi liquid theory requires be 
satisfied on the lowest energies \u\ < Ul, where lol is the 
low-energy scale characteristic of the coherent Fermi liquid 
state. The LMA itself is considered in §4, first in general 
terms applicable to an essentially arbitrary diagrammatic 
approximation for the inherent two-self-energies, and in- 
cluding the issue of symmetry restoration that is central 
to the approach. The specific non-perturbative approxi- 
mation to the LMA self-energies implemented here is then 
discussed, together with the practical method of solution 
such that the dictates of both symmetry restoration and 
the Luttingcr theorem are satisfied. 



Results arising are presented in §5, with a natural em- 
phasis on the strongly correlated Kondo lattice regime. 
An overview of dynamics on all energy scales is first given 
(§5.1), encompassing both the 'low'-energy behaviour char- 
acteristic of the renormalizcd heavy electron state as well 
as non-universal energies on the order of bare bandwidths 
or the interaction U . In addition to illustrating the broad 
roles of asymmetry (in both the conduction band and /- 
levels) , and of lattice type, comparison is also made on this 
'all scales' level both to results for single-particle dynam- 
ics of the AIM (in which only a single correlated /-level 
is coupled to the conduction band), and to dynamics of 
the PAM arising at the crude level of static mean-field. 
In §5.2 the material dependence of the low-energy lattice 
coherence scale lol on bare model parameters is obtained, 
in the strong coupling Kondo lattice regime where lol is 
exponentially small; and its behaviour compared in turn 
to corresponding LMA results for the AIM Kondo scale 
u>k- The central issues of scaling are considered in §5.3: 
the resultant universal scaling behaviour of dynamics in 
terms of u/ = uj/ul, on all lu' scales, and including the 
dependence of scaling dynamics on conduction band fill- 
ing. At low-u/ the scaling spectra exhibit coherent Fermi 
liquid behaviour, crossing over with increasing energy to 
logarithmically slow spectral tails. The latter are found to 
be independent of both conduction band filling and lat- 
tice type, and to have precisely the same scaling form as 
those for an AIM; establishing thereby the crossover from 
low-energy coherent Fermi liquid behaviour to effective 
incoherent single-impurity physics on high-o/ scales, but 
still in the w/a>£-scaling regime. A discussion of results 
obtained here in relation to the issue of two-scale 'exhaus- 
tion' physics [37,38] is given in §5.4; and some concluding 
remarks are made in §6. 
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2 Background 

The Hamiltonian for the PAM, H 
given in standard notation by: 



H r + Ht 



H, 



hybi 



H = zZ ec 4r c «7 - l zZ c « 

»,<r (*,.?),<' 



cr C jff 



+ E ( £ / + %fl-Ji-<r) flfia +Vj2(fl C Ur + ^) 

(2.1) 

The first two terms describe the uncorrelated conduction 
(c) band, H c ; with c-orbital site energies e c and nearest 
neighbour hoppings Uj = t, rescaled as t ex t*/\/Z c in the 
large dimensional limit where the coordination number 
Z c — > oo [3-6] (with t* the basic energy unit). The second 
term, Hf, describes the correlated /-levels, with site en- 
ergies ef and on-site repulsion Uff — U; while the final 
term, Hh y b couples the c- and /- subsystems via the local 
hybridization matrix clement V. Throughout the paper 
the Fermi level is taken as the origin of energy, u>p = 0. 

The model is thus characterized by four independent 
'bare'/material parameters, namely e c /i*, V/£*,e//i* and 
U/t* (with t* = 1 taken throughout) - a huge param- 
eter space in comparison e.g. to the Hubbard model. In 
previous LMA work on the PAM [34,35] we have stud- 
ied the particle-hole symmetric model appropriate to the 
Kondo insulating state; for which e/ = — -j and e c = 

0, with consequent occupancies nf = ^^ < fj a f icr > 

and n c = J2 a <cj <T c ilJ > = 1 for all U. Here we consider 
the generic asymmetric case, encompassing heavy Fermion 
metals (and with the symmetric PAM recovered as a par- 
ticular limit). Particle- hole asymmetry itself enters the 
problem in two ways [8] . (a) Conduction band asymme- 
try, reflected in e c ^ which, as detailed below, specifics 
the centre of gravity of the free (V = 0) conduction band 
relative to the Fermi level, (b) /-level asymmetry which, 
as for an impurity Anderson model [27], is embodied in 
the parameter 



1 



'/ 



1 



2f/ 

U 



(2.2) 



such that r\ = corresponds to particle-hole symmetric /- 
levels. The bare parameter set may thus be taken cquiva- 
lently as e c , V, U and r\. We shall find this choice to be con- 
venient in the following (and in fact necessary to describe 
universal scaling behaviour in the Kondo lattice regime, 
see §5.3). 

While the LMA developed here encompasses all inter- 
action strengths, the regime of primary physical interest 
is of course that of the strongly correlated Kondo lattice 
(KL): nf — ► 1, but with arbitrary conduction band fill- 
ing n c . The underlying low-energy model, obtained from 
the PAM to leading order in V 2 , is a Kondo lattice model 
(KLM); the KL regime of the PAM arising when e/ = 
-|e/| for |e/|/A) > 1 and (U - \e f \)/A > 1, where 
A = TrV 2 d^(0) (with 4(0) the free (V = 0) conduction 



band density of states at the Fermi level) . The approach 
to the KL is not therefore unique, in that ny — ► 1 arises 
for any asymmetry r/ = 1 — 2|e/|/J7 G [— 1, 1] on progres- 
sively increasing the interaction U. This is reflected in the 
fact that the associated KLM contains in general both ex- 
change and potential scattering contributions (the latter 
vanishes as the asymmetry 77 — > and is omitted in most 
studies of the KLM per se, which thus correspond to sym- 
metric /-levels alone but with asymmetry retained in the 
conduction band). 

Granted even a dominant interest is the strong cou- 
pling KL regime, the resultant bare parameter space of 
the PAM (or KLM) nonetheless remains 'large', as above. 
The KL regime is however characterized by a low-energy 
lattice scale, u>l, diminishing progressively with increas- 
ing interaction strength and expected to be exponentially 
small in strong coupling [7-25]. This scale is of course it- 
self a function of the bare material parameters; but that 
dependence is a subsidiary issue in comparison to the ex- 
pectation that physical properties of the PAM should ex- 
hibit universal scaling in terms of lo/lul and/or T/u>l, in a 
manner largely independent of the bare parameters them- 
selves. Understanding aspects of such scaling behaviour 
will be a central theme of the present work. 

Our specific focus in this paper is on local single-particle 
dynamics of the T — PAM, embodied in G u (uj) <-> 
G f u (t) = -i < f(f lrT (t)fl) > and likewise G?» for 
the c-levels, with corresponding local spectra D^(u>) = 
— isgn((jj)ImG^(o>) (and v = c or /). 

We begin with some remarks on the trivial limit V = 
where (Eq.(2.1)) the /-levels decouple from the free con- 
duction band. The latter is specified by its local propaga- 
tor denoted by go(<-o), with corresponding density of states 
(dos) dg(w) (= iV -1 X) Q <J(w - e a ) with e.g. e a = e k for 
a Bloch decomposible lattice); and it will prove useful to 
denote by H(z) the Hilbert transform 



/oo 
de 
-CO 



A)(e) 



(2.3) 



for arbitrary complex z, where po(u>) = <io(w;e c = 0) 
denotes the free conduction band dos for e c = (see 
Eq.(2.1)). The free c-electron propagator g§(w) is then 
given by 



5(w) = H(^+ - e c ) 



[u> + -e c - S (u>)] 



(2.4a) 
(2.4b) 



with uj + = lu + isgn(uj)0 + here and throughout, such that 
d^(uj) — pv{uj — e c ). Eq.(2.4b) simply defines the Feenberg 
self-energy [39,40] used below, with Sq(co) = S\go\ a func- 
tional of g§ alone (since g^ = H(5+l/gg)) from Eqs.(2.4)). 
While our subsequent discussion holds for an essentially 
arbitrary Po(lo) and hence band structure embodied in 
d()(co), specific results will later be given for the Bethe lat- 
tice (BL) and hypercubic lattice (HCL); for which within 
DMFT the normalized po (e) are respectively a semi-ellipse 
and an unbounded Gaussian, given explicitly (£* = 1) by 
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[3-6] 



Po(e) = -[l-e 2 ] 

A)(e) = ^exp(- 

V 71 " 



< 1 



BL 



HCL 



(2.5a) 
(2.5b) 



Since cIq(lo) — po(u> — e c ) is simply a rigid shift of po(w) 
(the free conduction band is non-interacting), conduction 
band asymmetry is thus embodied in e c itself, with e c = 
the symmetric limit. 

We turn now to the full local Green functions for the 
homogeneous paramagnetic phase of interest, for which 
the G^(lo) = G v (lo) are site-independent. The essential 
simplifying feature of DMFT - and the key aspect of it 
as an approximation to finite-dimensional systems - is 
that the /-electron self-energy is site-diagonal (momen- 
tum independent) [3-6] ; and from straightforward applica- 
tion of Feenberg renormalized perturbation theory [39,40] 
the G u (lo) are given by 

-l 



G°{u>) = 



G f (o 



V z 



LO+ — €f — Ef(uj) 

V 2 



-S(to) 



(2.6a) 



£/M 



CiT 



S(u>) 



[lo+ - e/ - Ef(w)] 



V 2 



[u>+ - e/ - E f (io] 



f G» 



(2.7) 

where Ef(co) = E?(u) — isgn(u>) EI(lo) is the conven- 
tional single self-energy (and the identity Eq.(2.7) follows 
from Eqs.(2.6)). The Feenberg self-energy S(lo) = S[G C ] is 
moreover precisely the same functional of the full G c (oj) as 
it is of 5 g(w) in the trivial limit of V = (e.g. S = \tlG c 
for the BL). In consequence, G c (lo) is given directly using 
Eqs.(2.6a,4,3) by 



G C H-H( 7 ) 



where 
7(w) 



V 2 



- e f - £>M ' 



(2.8) 



(2.9) 



For an arbitrary conduction band po(e), the approach 
to the full interacting problem is clear in principle: given 
the self-energy Ef(co), and hence 7(0;) from Eq.(2.9), G c (lo) 
= H(7) follows directly from Hilbert transformation and 
G' (lo) from Eq.(2.7). But practice is another matter: the 
hard part is to find a suitable approximation to the self- 
energy that can handle non-pcrturbativcly the strongly 
correlated physics of the KL regime, as well as the weak 
coupling regime of interactions (which itself is readily han- 
dled by plain perturbation theory or simple vari- 
ants thereof, see e.g. [14-17]). It is this impasse the LMA 
seeks to break, via use of an underlying two-self-energy 
description [26,27,34,35] as detailed in §4. In addition of 
course the problem must be solved iteratively and self- 
consistently, because an approximate Ef (lo) will itself be 
in general a functional of self-consistently determined prop- 
agators; that being a detail (albeit an important one) to 
which we likewise turn in §4.3,4. 



2.1 Non-interacting limit 

Before proceeding we comment briefly on the non- 
interacting (NI) limit, U = 0, the local propagators for 
which are denoted by Gq(lo) (or Gq(lo; e c ,ef,V 2 ) if ex- 
plicit dependence on the bare parameters is required) with 
corresponding spectra Dq(lo). Itself trivially soluble, the 
importance of the NI limit and rationale for discussing it, 
resides in its connection to the fully interacting problem 
via both Luttingcr's theorem [36] and the quasiparticle 
behaviour of the full D u (lo) at sufficiently low lo; as con- 
sidered in §3 below. The Gq(lo) are given by Eqs. (2.6-9) 
with Ef = and j(u>) — ► 7o(w) = lu + — e c — V 2 /[u> + — ef], 
with resultant spectra 



D§(w) = po(7o) =Po [u)-e 



V 2 



f(,.A — 



D> (u) 



V 2 



(we/) 2 



Po w - e c 



and hence total band filling n to t = n® + n f given by 



u)-e f 

V 2 
u-e f 

i° c + n° f 



(2.10a) 
(2.10b) 



(2 ' 6b) l(n° c + n° f ) =jl oo Mm») + D f (u)] 



L 



-ec+l/if 



p (e)de + 9(-~e f ) (2.11) 



where if = £f/V 2 (and 8(x) is the unit step function). 
For band fillings n to t & (0, 4) the system is generically 
metallic, with a non-zero Fermi level dos Dq(uj = 0). But 
since 70 diverges as lo — > e/, the spectral functions in the 
vicinity of lo — ef have the same behaviour as the tails of 
the bare conduction band po(e). So for a typical bounded 
Pa(e), e.g. the BL Eq.(2.5a), a hard spectral gap opens up 
in the neighbourhood of ey; while for an unbounded po(e), 
e.g. the HCL, a (strictly) soft gap arises at lo — ef. The 
system is of course insulating - a well known hybridization 
gap insulator [41] - only if the Fermi level (lo = 0) lies 
in the gap (excluding the trivial case of wholly empty or 
full bands); and from Eqs. (2. 10, 11) the sufficient condition 
for this to occur is readily seen to be ntot — 2, i.e. half- 
filling, which holds also for the fully interacting problem 
now considered. 



3 Adiabatic continuity, and quasiparticle 
behaviour. 

On increasing the interaction U from zero the system re- 
mains perturbatively connected to the NI limit; i.e. is a 
Fermi liquid, a statement applicable both to the metal- 
lic heavy Fermion (HF) state and the Kondo insulating 
(KI) phase which likewise evolves continuously from the 
non-interacting hybridization gap insulator [34,35]. This 
adiabatic continuity requires that the Luttinger integral 
vanish [2,36], i.e. that 



J-oc n ow 



(3.1) 
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What may be deduced on entirely general grounds 
from I i, = ? To that end note first that the local propa- 
gators G v (w) [y = c or /) may be expressed as 



G» 



dep (e)G»(e;u;). 



(3.2a) 



Thec-e\ectr<mG c (e;u)(=N- 1 Y /k [u + ~e c -e w -I] c (uj)}- 1 
with £k = e) follows directly from Eqs. (2.8,3) as 



G c (e;w) 



V 2 



e/ - S f (w) . 
while Gf(e;cu) follows in turn using Eq(2.7) as 

V 2 



Gf{e-Lu) 



fc / - ^1 



S f {w) 



uj^ - e c - e 



(3.2b) 



(3.2c) 



And in terms of the G v (e; u>) note that the total band 
filling n to t = n c + rif is given generally by 



/oo 
dep (e) 
-CO 



f ^[G c ( e ^) + G^( e ;c)] 



(3.3) 



Now use Eqs. (3. 2a, c) in Eq.(3.1), It = 0, together with 
the identity (from Eq.(3.2)) 



d^fju 



du> 



^-G f (e;u)= [G c (e;u) + G f (e;u)] 



du> 



In [(w + - e c - e)(w+ - e, - 27/ (w)) - F 2 



and perform the ^-integration. This yields 

f°° de 
i(n c + n/) = / — po(e)g(e) 



(3.4) 



using only that El(u> = 0) = 0, holding for both the HF 
and KI states; where 

g(e) = tan- 1 (r?(e + e c + e})/[e}(e + e c ) - V 2 ]) 

(with r/ = + ), and the renormalizcd level 



6} = e/ + £?(w = 0) 



(3.5) 



is thus defined. The e- integration in Eq.(3.4) is then read- 
ily performed to give the desired result 



\{n c + nf) 



-ec+l/?J 



p (e)de + 6{-? f ) (3.6) 



where e* f = eJ/V 2 . 

Eq.(3.6) is equivalently a statement of Luttinger's the- 
orem for the Fermi surface of the PAM, for the relevant 
case of a local, momentum independent self-energy appro- 
priate to DMFT (the Fermi surface is of course "large" , 
including /- and c- electrons, see also [38]). Three points 



should be noted about Eq.(3.6). (i) First and most impor- 
tantly we see it to be exact, following directly from II = 
without further approximation, (ii) It amounts to a simple 
renormalization of the NI limit result Eq.(2.11); being of 
just that form but with the bare level e/ = £f/V 2 replaced 
by the renormalized level e* f = [ej- + S^O^/V 2 ; which 
is thus determined via Eq.(3.6) for given filling ritot (and 
e c ). (iii) Eq.(3.6) is the direct analogue for the PAM of the 
Fricdcl sum rule for an AIM [2,42] , which relates the excess 
impurity charge (rn mp ) to the renormalized impurity level 
e* , and which likewise follows directly from 7^ = for 
the impurity model; see also §4.5. In physical terms that 
parallel is entirely natural, given the connection to an ef- 
fective impurity model which is inherent to DMFT [3-6] 
(Eq.(2.6b) for G' (u>) being of effective "single-impurity" 
form Gf(ui) = [uj + — ef — £f(u>) — Z\ e //(w)] _1 with an 
effective, w-dependent hybridization A e ff(ui) = V 2 [uj + — 
e c — S'((x')]- 1 ). Finally, we add that imposition of Eq.(3.6) 
as a self-consistency condition will play an important role 
in the LMA developed in §4ff. 

The second key implication of adiabatic continuity is 
that the limiting low-w behaviour of the propagators G v (u>) 
amount to a renormalization of the NI limit, which is of 
course the origin of the renormalized band picture [2,43]. 
This follows simply by employing the leading low-w ex- 
pansion of £f{u>), 



£/(w) 



Sf(0) 



1 

z- 1 



(3.7) 



withZ = [f — (<9i7?(u;)/da;) w= o] _1 the quasiparticle weight/ 
mass renormalization and SAui) neglected as ui — > (since 
£f((jj) oc to 2 for the HF metals or vanishes in the gap for 
the KI case). The low-w behaviour of the G v (uj) then fol- 
low from Eqs. (2.6) as 

G» ~Gg(w;e c ,Ze},ZV 2 ) = G c H (3.8a) 

G f (u>) ~ G f a (uj; e c , Ze* f , ZV 2 ) = Z&(u>) (3.8b) 

with Gq the NI propagators (§2.1) and the quasiparticle 
Green functions thus defined; with corresponding spectra 



D c {uj) ~ p (u) - e c - 



ZV 2 



(u - ZA) 



) = D c (co) 



D f (u) 



2x/2 



Z 2 V 

(c - Ze*) 2 P ° {uJ - ec 



ZV 2 



(u - Ze* 

zb f {uj) 



-) 



(3.9a) 



(3.9b) 



(e*f is the renormalized level, Eq.(3.5)). And the total band 
filling ntot = n c + nf, calculated from the quasiparticle 
propagators as \ntot — J_ dui [D c (lo) + D^ (ui)], correctly 
satisfies the exact result Eq.(3.6). 

Eqs(3.8,9) embody the quasiparticle behaviour of the 
PAM, and have important implications for the scaling be- 
haviour of D v (uj) in the strong coupling regime of pri- 
mary interest, as now considered. In the KL regime where 
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n/ — ► 1, the quasiparticle weight Z becomes exponen- 
tially small (as considered explicitly in §5.2). Defining a 
low-energy lattice coherence scale by (£* = 1) 



u L = ZV 2 



(3.10) 



the scaling behaviour of dynamics corresponds to con- 
sidering finite u/ = u)/u>l in the formal limit u>l — * 0. 
Eqs.(3.9) then yield 



connection between the KL regime of the PAM, and single- 
impurity Kondo physics, on suitably large energy and/or 
temperature scales. 

Before proceeding to the LMA we mention one fur- 
ther implication of Eq.(3.12) applicable to the KL regime: 
together with the exact result Eq.(3.6) it gives 



D c (us) ~ pa(-e c 



V 2 Df(u 



(u/ - e}) 2 



Po(-e c - 



) 



(u>> - ~e}) 



(3.11a) 



) (3.11b) 



where 'bare' factors of lo = lullu' may be neglected, and 
e*f = e* f /V 2 . Moreover, in the KL regime, e* f is solely de- 
pendent upon e c : from Eq.(3.11b) with ^n/ = J_ dioD$ (a. 
(Df( U ) = Df(u)/Z), 

f-0 



2 n / = 



do; 



1 



— OO 

6c + lA7 



(c'-?*)^o(- £ 



1 



depo(e) + #(-?;•) (3.12) 



whence e* = e*(e c ) as n/ — » 1 (and in addition sgn(ej) = 
sgn(e c )). 

Eqs.(3.11) embody the low-w behaviour of the single- 
particle spectra D"(u>), in the KL regime where nf — > 
1 but for arbitrary conduction band filing n c ; regarding 
which the following important points should be noted. 
(i) Eqs.(3.11) show that both D c (uj) and V 2 Di\io) (and 
not therefore Df(w) itself) exhibit one-parameter univer- 
sal scaling in terms of u>' = u/u>l- with no explicit de- 
pendence on the bare material parameters U, r\ (or e/) 
and V 2 ; and dependent solely upon e c (or equivalcntly 
on the conduction band filling n c = n c (e c ), see below) 
which itself determines the renormalizcd level e* f = e*Je c ) 
as above, (ii) Eqs.(3.11) provide explicitly the limiting be- 
haviour that, as \lu'\ — \u>\/u>l — * 0, must of necessity be 
recovered by any credible microscopic theory; and direct 
comparison of LMA results to which will be given in §5.3. 
Of equal importance however, the simple results above 
are asymptotically valid only as w' — ► 0, and prescribe 
neither the o/-range over which Eqs.(3.11) hold nor the 
general ^'-dependence of the scaling spectra - for which 
a real theory is required, (iii) The particle-hole symmetric 
PAM discussed in [34] (for which n/ = 1 = n c ) is just 
a particular case of the above, in which e c = and the 
renormalizcd level e* f = (by symmetry); and where the 
low-energy lattice scale ojl (Eq.(3.10)) is precisely the gap 
scale characteristic of the Kondo insulating state [34,35]. 
Finally, scaling arguments per se are obviously indepen- 
dent of how the low-energy KL scale ujl = Wi(e c , U, rj, V 2 ) 
itself depends upon the bare parameters; an issue of in- 
trinsic interest that has long attracted attention (see e.g. 
[8,21-25]) but which we believe (as argued in §5) to be 
in large part a red herring in understanding the expected 



depo(e) 



(3.13) 



for the c-band filling. This shows (a) that n c = n c (e c ) 
is indeed determined by e c as noted above; and (b) that 
the resultant n c is just that for the free (V — 0) conduc- 
tion band, for which (§2) rfg(w) = po( w — £c) with \n c = 

J_ dio d^{bj). In physical terms this is natural, since from 
Eq.(3.8a) the effective hybridization is ZV 2 = u>l, expo- 
nentially small in the KL regime such that the net con- 
duction band filling is in effect independent of coupling to 
the /-levels. We shall comment further on the latter in §5. 



4 Local Moment Approach 

The discussion thus far has been couched in terms of the 
single self-energy Sfiui) which, via diagrammatic pertur- 
bation theory in the interaction strength, provides the 
conventional route to dynamics. A determination of the 
propagators in this way is not however mandatory. Indeed 
while fine in principle there are good reasons to avoid it; 
stemming from the practical inability of conventional per- 
turbation theory, or partial resummations thereof, to han- 
dle the strongly correlated regime of primary interest. For 
this reason the LMA [26-35] takes a different route to the 
problem, employing instead a two-self-energy description 
that is a natural consequence of the mean-field description 
from which it starts. Here we first consider the implica- 
tions of such in general terms, independent of subsequent 
details of implementation (§4.3,4 ) and not confined to the 
symmetric PAM considered hitherto [34,35]. 

There are three essential elements to the LMA [26-35]. 
(i) Local moments ('/x'), regarded as the first effect of in- 
teractions, are introduced explicitly and self-consistently 
from the outset. The starting point is thus simple broken 
symmetry static mean- field (MF, i.e. unrestricted Hartree- 
Fock); containing two degenerate, locally symmetry bro- 
ken MF states corresponding to [i = ±|/u| [34]. While 
severely deficient by itself (see e.g. [26,27,31,34] and be- 
low), MF nonetheless provides a suitable starting point for 
a non-pcrturbativc many-body approach to the problem, 
(ii) To this end the LMA employs a two-self-energy de- 
scription that follows naturally from the underlying two 
local MF saddle points; with associated dynamical self- 
energies built diagrammatically from, and functionals of, 
the underlying MF propagators, (iii) The third and most 
important idea behind the LMA is that of symmetry restora- 
tion [27,28,32,34]: self-consistent restoration of the broken 
symmetry inherent at pure MF level; and in consequence, 
as discussed below, correct recovery of the local Fermi liq- 
uid behaviour that reflects adiabatic continuity in U to 
the non- interacting limit. 
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Within a two-self-energy description the local propa- 
gators G v {ui), which are correctly rotationally invariant, 
are expressed formally as (c/Eqs.(2.6)) 



where (with a = f/J, or +/—) 



G%(u) = [u> + -e c -- 



V 2 



-e/- E a {u) 



-S{uj)]- 1 (4.2a) 



Gt(u) = [u; + -ef-S a (u)- 



V 2 



ui+ — e c — S(u>) 



(4.2b) 



and S(u>) = S[G C ) as usual (the reader is referred to [34] 
for further, physically oriented discussion of these basic 
equations). The /-electron self-energies E a (u>) are conve- 
niently separated as 



E„(u>) = %{n - a\ji\) + E„{u) 



(4.3) 



where the first term represents the purely static Fock bub- 
ble diagram which alone is retained at pure MF level (with 
n and |/2| given explicitly by Eq.(4.12) below); and where 
the second term, E a {uj) — E^(u) — isgn(w)E I (u), is the 
key dynamical contribution mentioned above ('everything 
post-MF'). 

Eqs.(4.1,2) are the direct counterparts of the single 
self-energy equations Eqs. (2.6) (to which they would triv- 
ially reduce if E a {u>) = Ef(u>) for each a). For an arbitrary 
conduction band dos po(w) and any given {E a (u))}, they 
are likewise readily solved (cf the discussion of Eqs. (2.6a, 
9)): defining 



7o-(w) = u + - e c 



V 2 



(4.4) 



UJ+ - €f - E a {uj) 

such that G c {uj) = \Y, a b° - S ^ (Eqs. (4.1, 2a)), and 
comparing to G c (uj) = [7 — S]^ 1 (Eqs. (2. 6a, 9)), the 7 CT 's 
are related to the single 7(0;) (Eq.(2.9)) by 



7(w) = ±[7tM+7iM] + 



[1( 7t H-7|M)] 2 
S( w )-I[ 7T («)+ 7i ( w )]- 



(4.5) 



Given E a {u>) and hence 7 <t(w), this equation together with 



S(u) =7- 



H(i) 



(4.6) 



(from G c (lo) = [7 — S}^ 1 and Eq.(2.8)) may be solved 
straightforwardly and rapidly in an iterative fashion; em- 
ploying an initial 'startup' for S (typically S — \go(uS) 
with go the free conduction band propagator with dos 
po(u))). With S(w) then known, the G v (u>) follow directly 
from Eqs. (4. 1,2). 

The conventional single self-energy Ef(u>) follows im- 
mediately, essentially as a byproduct, because solution of 
Eqs. (4.5,6) determines both S(u) and 7(0;), whence (from 



Eq.(2.9)) E f (u>) = uj+ -e f -V 2 [uj+ -ec-jiuj)}- 1 follows; 
which relation may be expressed cquivalently as 



(4.1) ^/M = |[^ T (w)+A(w)]+ 



[i(£ T H-^H)] 2 

(4.7) 

Here Q(u>) is the usual host/medium propagator [44], given 
by Q-\w) = [{Gf (w))- 1 + Ef{u)\ - [w+ - e f - V 2 {uj+ - 
e c — S(lu)) x ] with corresponding spectral density T>(w) = 
— sgn(oj)Im^(a;) (and which in physical terms includes 
interactions on all sites other than the local site i {G{w) = 
Gu{uj)) [44]). The conventional single self-energy may thus 
be obtained directly given the {E a (u>)}, although obvi- 
ously not vice versa, and the underlying two-self-cnergy 
description may be viewed equivalently as a means to ob- 
tain Ef(uf). The particular class of diagrams retained in 
practice for the dynamical {E a (uj)} (see Eq.(4.3)) will be 
detailed in §4.3; at present none need be specified. 

At the pure MF level of unrestricted Hartee-Fock, dy- 
namical contributions to the E a (u)) are of course neglected 



entirely and E a (u)) = E a 



-(n - a\fj,\) (with the MF lo- 



cal /-level charge n = n and moment \\x\ = \\x\ determined 
in the usual simple fashion, §4.2). From Eq.(4.7) the single 
self-energy at MF level is then 



Ef F (Lu) = ^Un +7 ^ 



au\v\) 2 



Gr(w)-\Un 



(4.8) 



(with Gq{uj) the corresponding MF medium propagator, 
whose Fermi level spectral density T>q(uj = 0) is read- 
ily shown to be non-zero). From this the basic deficiency 
of pure MF is clear: if the local moment \fi\ ^ 0, then 
from Eq.(4.8) the Fermi level lm.Ef F {oj = 0) ^ and 
Fermi liquid behaviour is violated — wholly wrong, albeit 
arising naturally because the resultant degenerate MF lo- 
cal moment state is not perturbatively connected to the 
non-interacting limit. While this problem would not oc- 
cur if I /1 1 = were enforced a priori (restricted Hartree- 
Fock), another one then arises at post-MF level; for from 
Eq.(4.8) the two- and single- self-energy descriptions then 
coincide, with E¥ f (uj) = \Un merely the static Hartree 
contribution, producing a trivial energy shift to the non- 
interacting propagators. Subsequent construction of the 
dynamical Ef(u>) via conventional perturbation theory in 
U employing these propagators, is equivalent to expand- 
ing about the restricted Hartree-Fock saddle-point. But 
when local moments can form at MF level this singlc- 
determinantal saddle point, unlike its unrestricted MF 
counterpart, is unstable to particle-hole excitations. It is 
this in turn that is readily shown to underlie the familiar 
divergences arising within conventional perturbation the- 
ory if one attempts to perform the 'natural' diagrammatic 
resummations (such as RPA) that one expects physically 
are required to capture regimes of strong electronic cor- 
relation, and the general inability to surmount which has 
been a plague on all our houses [2,45]. 

The LMA seeks to surmount these problems by (a) 
retaining the two-self-energy description, with the inher- 
ent notion of local moments and essential stability of the 
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underlying MF state; while (b) incorporating many-body 
dynamics into the associated self-energies {^(w)} in a 
simple and tractable fashion, and in such a way that Fermi 
liquid behaviour is recovered at low-energies. 

4.1 Symmetry restoration 

This brings us to the key notion of symmetry restora- 
tion (SR), now sketched briefly in the generic context 
of heavy Fermion (HF) metals in the asymmetric PAM, 
where it arises from the obvious question: under what 
conditions on the {E a (uj)} will the /-electron single self- 
energy Sf(uj) exhibit Fermi liquid behaviour as w — > 0, i.e. 
will El(w) ~ 0(o> 2 )? This may be answered simply by em- 
ploying a general low-frequency Taylor expansion for the 
E a (uj) in Eq.(4.7), along precisely the same lines as in [27] 
for the Anderson impurity model. That is merely a matter 
of algebra, and from it one finds the necessary/sufficient 
condition for El(uj) ~ O(oj 2 ) is that 



at this stage to comment briefly on MF itself; denoting 
the MF propagators by g v a {w) (and g v {u) = ± £ ff #(w)). 
These follow from Eqs.(4.2) as 



£f(w = 0) = £f(w = 0). 



(4.9) 



Moreover, with Eq.(4.9) satisfied then from Eq.(4.7) (i) all 
self-energies coincide at the Fermi level, i.e. 

Ef (u> = 0) = E*(u> = 0) (4.10) 

for either spin a; (ii) the leading \ow-ut behaviour of E¥(u>) 
is as in Eq.(3.7), with the quasiparticle weight Z = [1 — 
{dEf(w)/du)) IJJ= o)}~ 1 given by Z- 1 = \{Z^+Z^) where 
Z„ = [1 - (<9£*(u;)/dw)u;=o] _1 is thus defined; and (iii) 
the quasiparticle behaviour embodied in Eqs.(3.8) for the 
propagators G u (u) is thus recovered. 

Eq.(4.9), a condition upon the E„(uj = 0) solely at the 
Fermi level, is the SR condition that is central to the LMA. 
It is quite general, being precisely the condition found 
for Anderson impurity models, whether metallic [26,27] or 
pseudogap AIMs [31,32]; and likewise for the particle-hole 
symmetric limit of the PAM [34] where it guarantees per- 
sistence of the insulating gap with increasing interaction 
strength, reflecting the 'insulating Fermi liquid' nature of 
the Kondo insulating state. The general consequences of 
SR are correspondingly common to all these problems: 
In practice, Eq.(4.9) amounts to a self-consistency equa- 
tion for the local moment |/x| (supplanting the pure MF 
condition for \/i\), see §4.3,4. Most importantly, §4.4 and 
[26-28,34], imposition of SR as a self-consistency condi- 
tion generates a non-vanishing low-energy spin-flip scale 
oj m , manifest in particular in the transverse spin polar- 
ization propagator, whose physical significance is that it 
sets a non- vanishing timescale, r ~ h/u) m , for the restora- 
tion of the broken spin-symmetry endemic to the pure 
MF level of description; and which in the present context 
is cquivalently the low-energy Kondo lattice scale, with 
w m (xw L = ZV 2 (Eq.(3.10)). 

4.2 Mean-field 

Since the self-energies E a (uj) are built diagrammatically 
from the underlying MF propagators, it is appropriate 






V 2 



e f 



uj+ — e/ + ax 

V 2 
- ax 



S^)}- 1 (4.11a) 
1 ~ 1 (4.11b) 



lu+ — e c — S(io) 

where S(u>) = S^ ]; and we have written the purely static 
e f + E a {w) = e f + E° as e f + E° = e f -ax, with x = \U\n\ 
and e/ given at pure MF level by e/ = e/ + iUn. For any 
given e/ and x, explicit solution of Eqs.(4.11) for the MF 
propagators g^(u>) = g%(u>;ef,x) follows directly in one 
shot as described above (Eqs.(4.4-6)). And at pure MF 
level, the local moment |/x| and charge n are found from 
the usual MF self-consistency conditions |/i| = |/2| and 
n = n; where \p,\ = \p,(ef,x)\ and n = n(ef,x) arc given 
generally by 

\p,\ = / duj [dl(cj;ef,x) — dAu>;ef,x)\ (4.12a) 

J — oo 

duj [df(u>;e f ,x)+d i i (oj;ef,x)} (4.12b) 



(such that the static Fock bubble diagram, appearing in 
Eq.(4.3) for E„(u), is given generally by ^-(n — a\p,\)). 
In practical terms here it is obviously most efficient to 
work with fixed e/ and x: from Eq.(4.12a), \fi\ = \p,\ yields 
immediately U = 2x/\)i\, whence with n — n Eq.(4.12b) 
likewise gives directly e/ = e/ — ¥n. One can of course 
choose equivalently to specify the bare parameters U and 
€f (or rj = 1 + 2ef/U) from the beginning — which simply 
requires iterative cycling of the pure MF self-consistency 
equations — whence Eqs.(4.12) determine the pure MF 
values for x = \U\^i\ and e/. Results arising at pure MF 
level will be shown explicitly in §5.1. 



4.3 LMA: practice 

Beyond the crude level of pure MF it is of course the dy- 
namical contributions to the self-energies, E a (uS) (Eq.(4.3)), 
that are all important; and since the E a {bj) arc functionals 
of the underlying MF /-electron propagators {gl{u; e/, x)} 
themselves given by Eq.(4.11b), the E a (u)) = E„(w; ef,x) 
thus depend upon e/ and x. Hence, independently of the 
particular class of diagrams retained in practice for the 
dynamical self-energies, the question arises as to how x = 
■%U\n\ and e/ are determined in general (for any given set 
of bare model parameters)? To do so clearly requires two 
conditions. As discussed in §4.1, the symmetry restoration 
condition Eq.(4.9) must of necessity be satisfied in a Fermi 
liquid phase; and using Eq.(4.3) it may be cast as 

Ef(u = 0;e f ,x)-Ef(u) = 0;ef,x) =U\p,(e f ,x)\. (4.13) 

Likewise, as discussed in §3, adiabatic continuity requires 
that the Luttinger integral theorem Eq.(3.1) be satisfied, 



N.S.Vidhyadhiraja and David. E.Logan: Dynamics and scaling in the periodic Anderson model 



i.e. 



I L (ef,x) = Im 



du> dSf(u>) 



du> 



G f {u 







(4.14) 



where the Luttinger integral itself depends necessarily on 
e/ and x. The essential point is obvious: these two equa- 
tions are sufficient to determine x = 2^1mI an d e f m 
the general case, and in effect supplant the correspond- 
ing pure MF self-consistency conditions discussed above 
(recall that since the broken symmetry MF state itself is 
not adiabatically connected to the non-interacting limit, 
MF fails in general to satisfy either symmetry restoration 
or the Luttinger theorem). The optimal method for solv- 
ing these equations will naturally depend on the particular 
approximation employed for the dynamical {E a (uj)}; but 
that is an algorithmic detail, to which we return below. 
One further comment should be added here. It is readily 
shown that the particle-hole symmetric PAM considered 
in [34,35] (for which e/ = — -j,e c = 0, see §2), corre- 
sponds necessarily to e/ = 0, and that the Luttinger the- 
orem Eq.(4.14) is automatically satisfied by particle-hole 
symmetry. In that case solely the symmetry restoration 
condition is therefore required to determine x = t^ImI 
and hence the local moment, |/i|; precisely as employed in 
previous LMA work on Kondo insulators, [34,35]. 

While the preceding discussion is general, the final task 
is to specify the class of diagrams contributing to the dy- 
namical /-electron self-energies {E a (uj)} that we here re- 
tain in practice. These have the same functional form cm- 
ployed in [34,35] for the symmetric PAM, and may be 
cast as shown in Fig.l. The wavy line denotes the local 
interaction U , the double line denotes the broken sym- 
metry host/medium propagator Q-o-(ui) specified below 
(Eq.(4.16)); and the local /-level transverse spin polar- 
ization propagator 77~ CTCT (o>), likewise specified below, is 
shown as hatched. The diagrams translate to 

/OO 7 

=lg-„(u-ui)n-'"'(ui) (4.i5) 
-oo 27T1 

and retention of them is motivated on physical grounds, 
for they describe correlated spin-flip processes that are es- 
sential in particular to capture the strong coupling Kondo 
lattice regime: in which having, say, added a cr-spin elec- 
tron to a —cr-spin occupied /-level on lattice site i, the 
—cr-spin hops off the /-level generating an on-site spin- 
flip (with dynamics reflected in the polarization propa- 
gator n~ ac7 (oj)). The —cr-spin electron then propagates 
through the lattice in a correlated fashion, interacting fully 
with /-electrons on any site j =/= i (as embodied in the 
host/medium Q_ a (yj)); before returning at a later time to 
the original site i, whence the originally added cr-spin is 
removed (and which process simultaneously restores the 
spin- flip on site i). 

The renormalized /-electron medium propagator 
Q-„(ui), which embodies correlated propagation of an /- 
electron through the lattice, is given explicitly by (c/ its 
counterpart Q(u)) arising in Eq.(4.7)) 



G-a(u) = [^ + - e/ 



V 2 



y-rj\s\rKj\s\r\r\ra 



£ a («) 



= -G,\ 



O/i 



I'-C 



Fig. 1. Diagrams retained in practice for the dynamical /- 
electron self-energies £ a (u)). Wavy line: interaction U. Dou- 
ble line: renormalized host/medium propagator (see text). The 
transverse spin polarization propagator is shown hatched. 



with corresponding spectral density T>- a {w). Physically, 
Q-a-(u)) = Gu;-(t(w) includes interactions on all sites other 
than i (on which interactions occur at MF level); and the 
dependence of S a (uj) (Fig.l) on which accounts in effect 
for the hard-core boson nature of the spin-flips [34,46]. 
Diagrammatic expansion of Q- a {ui) in terms of MF prop- 
agators and self-energy insertions E^ a (w), and hence the 
infinite set of diagrams implicit in Fig.l for S a (uj), is dis- 
cussed further in [34,46] to which the reader is referred. 

The local (site-diagonal) polarization propagator en- 
tering Eq.(4.15) for U a (u>) is given at its simplest level 
by an RPA-likc particle-hole ladder sum in the transverse 
spin channel, namely 



n- 



r M 



u n~ 



r H[i-[/°7T- CT »r 



(4.17) 



with the corresponding bare polarization bubble °77 (u>) 
= °77 (u>; e/, x) expressed in terms of the broken sym- 
metry MF propagators {gf(u>; e/,x)}; referred to in [34] 
as 'LMAF. [Alternatively, one may readily renormalize the 
polarization bubbles in terms of the host/medium prop- 
agators {G„(w)}, so-called 'LMAIF [34]. Results arising 
from the two are however very similar [34], so we largely 
confine our attention in the present paper to LMA I, ex- 
cepting the explicit comparison between the two made 
in Fig. 9 below.] The °77 aa (uj) and hence II~ aa (uj)) are 
moreover readily shown to be related by II~ c " t (uj) = 
/7 <J ~ <J (— lo) [27,34]; whence only one such need be consid- 
ered explicitly, say 77 H (u>). Using this, and the Hilbert 
transform for 77 -1 (u), Eq.(4.15) for the dynamical self- 
energy reduces to 

/■OO J 

Z a (cj) = U 2 — Im/7+-(wi) fl^wijai^w + crwi) 



S(a 



- 1 (4-16) d + 



+ 9{-(JUJ l )g+ a {uJ + (TU Jl ) (4.18) 

where G ± a (oj) = J^ doo 1 f>_ a (u) 1 )6(±u) 1 )[u) - wi ± W+]" 1 
denote the one-sided Hilbert transforms such that Q_ a {uj) = 



10 
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4.4 Solution 

We now summarise the preceding discussion from the view- 
point of practical solution, and specify what we find to be 
a numerically efficient algorithm to solve the basic LMA- 
DMFT equations. 

The self-energies {£ a (u>)} are given in their entirety by 
Eq.(4.3), with the static Fock contributions |/S(e/, x)\ and 
n(ef,x) from Eq.(4.12). The dynamical contribution to 
the self-energy, S a (oo), is given by Eq.(4.18), with the po- 
larization propagator therein specified by Eq.(4.17); and 
the host/medium propagator Q_ a (ui) given by Eq.(4.16) 
(itself dependent on the Feenberg self-energy S{u) = S[G C ] , 
requiring as such an iterative, self-consistent solution of 
the problem). For given S a {ui), Eqs. (4.4-6) are the key 
equations to solve (as there discussed) for G c (uj) and S(ui); 
G* (uj) then follows directly from Eqs. (4. 1,2b). In addi- 
tion, centrally, both the symmetry restoration condition 
for S„{u = 0) (Eq.(4.10) or (4.13)) and the Luttinger 
integral theorem Eq.(4.14) — or equivalently Eq.(3.6) - 
must also be satisfied; which conditions, for given bare pa- 
rameters {e c ,V 2 ,ef,U}, determine both x = ^U\fi\ (and 
hence the local moment |/x|) and e/ that prescribe the un- 
derlying MF propagators. In this regard we note for use 
below that Eq.(4.4) for 7o-(w) may be written equivalently 
as 



7<tM 



V 2 



- e* - \E a {w) - 27,(0)] 



(4.19) 



where el (Eq.(3.5)) is the renormalized level, el — e/ + 
Sf(0) = ej + 27^(0) (for either spin a, as follows directly 
from symmetry restoration Eq.(4.10)); and from Eq.(4.3) 
we have used trivially that E a (ui) — E a (0) = E a (ui) — 

The particular algorithm employed is now specified, for 
an arbitrary conduction band po(e). In practice we choose 
to work with specified e c ,V 2 ,x and e/, with the bare pa- 
rameters U and €f determined by solution; rather than 
with e c , V 2 , e/, U specified and x, e/ then determined. The 
two are of course entirely equivalent; we simply find the 
former to be optimal in practice. So for any given e c , V 2 , e/ 
and x — |C/|/i|, the algorithm is as follows: 

(i) 'Startup'. Eqs. (4. 11) are first solved for the MF 
propagators {g%.(u>; e/, x)} (v = c or /), following the 
procedure specified in Eqs. (4.4-6). From this, the polar- 
ization bubble °77 (uj) (and hence 77 H (uj)) follows di- 
rectly, see Eq.(4.17). E a (ui) is given by Eq.(4.18), in which 
the host/medium propagator Q_ a (uj) is initially taken to 

be the MF gL a (uj), thus generating the 'startup' E a (ui). 

(ii) Symmetry restoration. The u> = SR condition 
Eq.(4.13) is now solved for the interaction U. This sim- 
ply requires varying U in Eq.(4.18) for E a (uj = 0) until 
Eq.(4.13) is satisfied (the [/-dependence of 27^(0) arising 
both from the trivial U 2 prefactor in Eq.(4.18) and the ex- 
plicit [/-dependence of 77 H (ui), see Eq.(4.17)). The local 
moment follows immediately, \fi\ = 2x/U. 

(iii) Luttinger condition. With an input guess for the 
renormalized level el, and hence 7 CT (w) from Eq.(4.19), 



Eqs. (4.5-6) are readily solved (as there described) for G c (u>) 
and S(co); and G* (uj) follows directly from Eqs. (4. 1,2b). 
The total band filling is trivially computed from h(n c + 

rif) = J_ dui[D c (u>) + D* (u>)\, and compared to the Lut- 
tinger condition Eq.(3.6) (in which e* f — e* f /V 2 ). The 
renormalized level el is then simply varied until Eq.(3.6) 
is self-consistently satisfied; the corresponding bare e/ fol- 
lows directly from e/ = el — 27^(0). 

(iv) The resultant S(u>) is then used in Eq.(4.16) to 
generate a new host/medium propagator Q- a (ui) ; and hence 
from Eq.(4.18) a new S a {uf). Now return to step (ii) and 
iterate until full self-consistency is reached. 

We find the above algorithm to be efficient, converg- 
ing typically after ~ 6 iterations and computationally fast 
on a modest PC. The outcome is a fully self-consistent 
solution of the problem, with {e c , V 2 , e/, U} and e/, x all 
known (uniquely so in practice) . If one wishes instead to 
work with a specified bare parameter set {e c ,V 2 ,ef,U} 
one simply repeats the above procedure, varying e/ and 
x = 2 C^ImI until the desired e/, U are obtained. The particle- 
hole symmetric PAM studied in [34,35], with e/ = — ^ and 
e c = 0, is a special case of the above algorithm; here e/ = 
(and likewise e% — 0), and the Luttinger condition is au- 
tomatically satisfied so that step (iii) above is redundant. 
Results arising from the fully self-consistent solution will 
be discussed in the following sections. 

Before proceeding we comment on the low-w behaviour 
of the transverse spin polarization propagator 77 ^ (uj), 
given by Eq.(4.17) and entering the self-energy Eq.(4.18). 
As mentioned in §4.1, this is characterised by a low-energy 
spin-flip scale denoted by ui m (and defined conveniently 
as the location of the maximum in ImT?" 1 {ui)). Such be- 
haviour arises in all problems studied thus far within the 
LMA [26-35] and has a common origin now briefly ex- 
plained. If the local moment |/z| had its pure MF value 
- i.e. if \fi\ was determined from the usual MF self- 
consistency condition (§4.2) \/j,\ = \p,{ef,x = \U\n\)\ with 
|/2(e/,x)| given generally by Eq. (4. 12a) — then it is straight- 
forwardly shown (following e.g. [27,34]) that 77 H (ui) given 
by Eq.(4.17) contains a pole at uj = identically. In phys- 
ical terms this reflects simply the fact that the pure MF 
state is, locally, a symmetry broken degenerate doublet, 
with zero energy cost to flip an /-level spin. This is cor- 
rect only in the 'free lattice' limit of vanishing hybridiza- 
tion V where the /-levels decouple from the conduction 
band, resulting in a degenerate local moment state (a 
limit that we note is recovered exactly by the LMA, non- 
trivially so from the perspective of conventional pertur- 
bative approaches to the PAM). Such 'cost free' spin- flip 
physics is not of course correct for the Fermi liquid phase 
of the PAM that is adiabatically connected to the non- 
interacting limit. But neither does it occur in this case, for 
the existence of an u> = spin-flip pole is readily shown to 
be specific solely to the pure MF level of self-consistency 
(i.e. arises only if |/z| has its pure MF value specified 
above). The key point is that within the LMA the lo- 
cal moment |/i| is determined from the symmetry restora- 
tion condition Eq.(4.10) (as in step (ii) above), which it- 
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self reflects adiabatic continuity (see §4.1). In consequence 
Imi7^ (u>) contains not an u) = spin-flip pole, but rather 
a low-energy resonance centred on a non-zero frequency 
u) m , whose physical content in setting the timescalc for 
symmetry restoration has already been noted in §4.1 (and 
which is equivalently the low-energy Kondo lattice scale 
cj m ocuj L =ZV 2 (Eq.(3.10))). 



4.5 Single-impurity model 

For sufficiently low energies and/or temperatures the be- 
haviour of the PAM is that of a coherent Fermi liquid [1,2], 
reflecting the lattice periodicity and embodied in the lat- 
tice coherence scale lol- However with increasing energy 
and/or temperature, it has of course long been known that 
a crossover should occur to an incoherent regime of ef- 
fective single- impurity physics [1,2]. For that reason it is 
traditional in studies of the PAM/KLM to compare to cor- 
responding results for the Anderson impurity (or Kondo) 
model; in which only a single /-level is coupled to the 
host conduction band, but otherwise with the same bare 
parameters as the PAM itself, (e c , V 2 ,Ef,U). Here we com- 
ment briefly on the LMA for the relevant Anderson im- 
purity model (AIM) itself [26,27], which will be used in 

§5- 

The local impurity Green function for the AIM, which 
we continue to denote by G*(ui), is given by G' (u>) = 



■£„G£(w) ( c/ Eq.(4.1)) where 



Gfa) 



e f -S a {u)-A{u)]-\ 



(4.20) 



One-electron coupling between the impurity /-level and 
the host is as usual embodied in the hybridization func- 
tion A(u>) = Ar(lo) — isgn(u;)Ar(u;), and is given sim- 
ply by A(u>) = V 2 g$(io); where g%{(J) is the free lattice 
(V = 0) conduction electron propagator (Eq.(2.4)) with 
corresponding dos rfg(w) = pa(uj — e c ) (and po{oS) given e.g. 
for the hypercubic and Bcthc lattices by Eqs.(2.5)). Note 
then that the AIM hybridization is both w-dependent and, 
for generic e c 7^ 0, asymmetric in lo about the Fermi level 
lo = (in contrast e.g. to the commonly considered wide 
flat-band AIM [2] for which A R (ui) = and Ar(w) = con- 
stant). This w-dependence in A{lo) will of course be ap- 
parent in AIM single-particle dynamics on non-universal 
energy scales (as seen e.g. in Fig. 5 below). But in the 
strong coupling Kondo regime of the AIM, characterised 
by an exponentially small Kondo scale lok — > 0, universal 
scaling behaviour of dynamics arises in terms oIlo /lok (see 
e.g. [26-28,33]). In this regime the w-dependence of A(u) is 
naturally immaterial, and only A(u> — 0) is relevant. With 
that in mind, for later use we denote by Aq — Ai{lo = 0) 
the local hybridization strength at the Fermi level, 



A = irV 2 p (-e c ). 



(4.21) 



The self-energies {E a (w)} for the AIM in Eq.(4.20) are 
again given by Eq.(4.3); with n and |/2| by Eq.(4.12), where 
now the MF spectral densities dl(u>) = d* a (uj;ef,x) nat- 
urally pertain to the MF propagators for the AIM, given 



by (cf Eq.(4.11b)) gf{u) = [u+ - e, + ax- A{lo)]-\ 
The dynamical contributions to the self-energies, {£&(lo)} 
(Eq.(4.3)) arc likewise given [26,27] by Eqs.(4.15) or (4.18), 
with the self-consistent host/medium propagator Q_ a {uj) 
appropriate to the PAM now replaced simply by the AIM 
propagator g_ a (ui). And the symmetry restoration condi- 
tion for the AIM is again given by Eq.(4.9) [27], whence 
£f(0) — -^(O) (for either spin a) where Sj(lo) here 
denotes the AIM single self-energy Finally, with e* = 

e/+I7?(0) denoting the impurity renormalized level (as for 
the PAM, Eq.(3.5)), the Luttinger integral theorem i£=0 
(Eq.(3.1)) yields directly the Freidel sum rule for the AIM 
[2,42]: e* mp + A R (0) = Z\ tan(f n imp ), with n imp as usual 
the excess charge induced by addition of the impurity, and 
which is the AIM analogue of Eq.(3.6) appropriate to the 
PAM. The LMA for the single-impurity model is read- 
ily implemented, as detailed in [27], with both symmetry 
restoration and the Luttinger theorem satisfied. 



5 Results 

We turn now to results arising from the LMA specified 
above. Following consideration of dynamics on all energy 
scales (§5.1), the dependence of the coherence scale lol 
on bare/material parameters is obtained in §5.2 and com- 
pared to corresponding results for the Kondo scale lok 
of the single-impurity Anderson model. The central is- 
sues are considered in §5.3: the u>/aj£-scaling behaviour of 
single-particle spectra in the strong coupling Kondo lattice 
regime, and their evolution from the low-energy physics 
characteristic of the coherent Fermi liquid through to the 
emergence at high energies of single-impurity scaling be- 
haviour. Finally, §5.4 discusses our results in the context of 
Nozieres' problem of 'exhaustion' [37,38] and recent work 
on that issue. 



5.1 All scales overview 

For obvious physical reasons the primary interest in the 
PAM resides both in the strongly correlated Kondo lattice 
regime, and on energies on the order of the coherence scale 
u>l = ZV 2 and (essentially arbitrary) multiples thereof. 
We begin however with an overview on all energy scales 
— encompassing 'band scales' lo ~ O(l) (1 = £*) and en- 
ergies lo r~j e/ or €f + U characteristic of the /-electron 
Hubbard satellites. In contrast to the low-energy sector, 
dynamics here will naturally be non-universal: dependent 
on essentially all bare material/model parameters, and lat- 
tice specific. An overview is nonetheless instructive, show- 
ing clearly the roles of asymmetry (in both the conduction 
band and /-levels), and of the lattice type, as well as qual- 
itative effects of depleting the conduction band filling. In 
addition, it enables broad comparison both to dynamics 
arising at the crude level of pure MF (§4.2) and to corre- 
sponding results for the Anderson impurity model (§4.5). 
Figs. 2 and 3 show spectra typical of metallic heavy 
fermion behaviour in strong coupling: U ~ 5.1 (x = ^?7|£i| 
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Fig. 2. All scales view of LMA c- and /-electron spectra (solid lines): D c (uj) and D* (u>) i 
panels) and hypercubic lattice (right panels). For U ~ 5.1 (x — |f/|/i| = 2.5), V 2 = 0.2, 
Corresponding MF spectra are also shown (dashed lines). 



U (= u)/t- M ), for Bethe lattice (left 
c = 0.3 and n = 1 - 2\e f \/U = 0. 



= 2.5), V 2 = 0.2 and e c = 0.3. In Fig.2, rj = l + 2e f /U = 
is taken — corresponding to symmetric /-levels e/ = — If j 
but with asymmetry in the conduction band (e c ^ 0). To 
illustrate the effects of the lattice, the c- and /-electron 
spectra are each shown for both the hypercubic lattice 
(HCL) and Bethe lattice (BL). In either case the /-level 
charge nj ~ 0.99. The conduction band fillings, likewise 
determined from spectral integration, differ little for the 
two lattices, n c ~ 0.64 (BL) and 0.69 (HCL) (each be- 
ing within ~ 2% of the asymptotic strong coupling result 
Eq.(3.13) for n c ). 

The overwhelming intensity of the /-spectra shown in 
Fig.2 is naturally in the Hubbard satellites, well sepa- 
rated from the band scales u> ~ 0(1) (and in consequence 
sharply distributed). Their peak maxima are symmetri- 
cally positioned about the Fermi level — reflecting the 
fact that the /-levels themselves are symmetric (rj = 0) 
- and largely unaffected by the presence of asymmetry 
in the conduction band. The most important feature of 
the /-spectra is of course the well known many-body res- 
onance at low energies. Its rich structure, considered in 
detail in §5.3, is naturally not resolved here. What is how- 
ever evident from Fig.2 is the relative unimportance of 
the host lattice in determining the /-spectra on the all 
scales level shown. This is in contrast to the local conduc- 



tion electron spectra (top panels, Fig.2). Here, aside from 
weakly remnant Hubbard satellites whose intensity dimin- 
ishes steadily with increasing U, the c-spectra are clearly 
dominated by the asymmetrically distributed envelope of 
the free (V = 0) conduction band spectrum, semi-elliptic 
for the BL and Gaussian for the HCL. As mentioned at 
the end of §3 this is physically natural, reflecting that in 
the strongly correlated regime the conduction band is very 
weakly coupled to the /-levels; albeit that such coupling 
is of course the key feature of the problem on low energy 
scales, where it leads to many-body structure in the D c (tu) 
(again barely visible on the scales shown). In Fig. 3, shown 
for the HCL with 77 = 0.4, there is now particle-hole asym- 
metry in the /-levels as well as in the conduction band; 
producing the additional spectral signature of asymmetry 
in the positions of the Hubbard satellites, but otherwise 
little change in broad terms. 

Figs. 2 and 3 also show direct comparison to the corre- 
sponding spectra at pure mean- field level (dashed lines). 
At first sight, and on the all scales level shown, these ap- 
pear to provide a reasonable first approximation to dy- 
namics. That is not of course the case on the all important 
low-energy scales that dominate the physics of the PAM in 
strong coupling: MF clearly lacks any hint of the many- 
body resonance in the /-spectra and its counterpart in 
D c (to) — unsurprisingly given the absence of correlated 
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Fig. 3. Local c- and /-electron spectra for the HCL: D c (cj) and D' (uj) vs u), shown for r\ 
parameters as in Fig. 2. Dashed lines: corresponding MF spectra. 



2| e/|/[/ = 0.4 with remaining 



electron dynamics at this crude level — and in fact for 
Df(uj) is seen to be qualitatively deficient for essentially 
all \lu\ < 1. For the c-clcctron spectra however, and again 
excepting the lowest energies, MF is qualitatively reason- 
able, reducing in strong coupling to precisely the free con- 
duction band spectrum (as is readily shown directly from 
Eqs.(4.11)). In addition, MF also captures qualitatively 
the dominant Hubbard satellites in the /-electron spectra; 
albeit that many-body broadening effects, arising from the 
spin-flip dynamics included in the LMA, lead both to a 
broadening and slight shift of the satellites (which can be 
understood quantitatively in terms of the w-dependence 
of the dynamical self-energies E a (ui), although we do not 
pursue that further here). 

We consider now the qualitative effect of depleting the 
conduction band filling n c , obtained by increasing e c sig- 
nificantly. Fig. 4 shows HCL spectra for e c = 1.5, U ~ 7 
(x = 3.5), f] = 0.8 and V 2 = 1.25; for which the resultant 
n c ~ 0.2 (and n/ ~ 0.85). Save for the large e c the re- 
maining parameters have no special significance, and the 
large V 2 has simply been chosen so that the resultant low- 
energy scale wi (discussed in detail in §5.2) is not so small 
as to be in effect invisible in the figure shown. Depleting n c 
in this way has a marked effect on the conduction electron 
spectra. In contrast to Figs. 2, 3 for e c = 0.3 — where the 
low-energy 'antiresonance' in D c (uo) is carved out of the 
free conduction band envelope — Fig. 4 shows that the 
low-n c conduction spectrum now contains a sharp low- 
energy resonance lying in the tail of the free conduction 
band envelope, akin to that appearing ubiquitously in the 
local /-spectra. The essential origin of the resonance is 
readily seen from the quasiparticle behaviour discussed in 
§3: from the quasiparticle D c (ui) Eq.(3.11a), the Fermi 
level D c {lo = 0) ~ Po{-e c + 1/e}) (with e* f = e* f /V 2 
and e*f the renormalized level); while for large e c = |e c |, 
Eq.(3.12) shows that — e c + l/e* f — ► as n/ — > 1. In con- 
sequence, D c (uj = 0) ~ po(uj = 0) is effectively pinned at 
the Fermi level to its free lattice limit (with po(0) = l/vA" 
for the HCL marked explicitly on Fig. 4); and the width 



of the resultant resonance is 0(ujl), as again follows from 
the quasiparticle D c (lu), Eq.(3.11a). The low-energy reso- 
nance arising for low n c in the c-electron spectrum thus re- 
flects directly the Fermi liquid nature of the ground state. 




Fig. 4. Depleting the conduction band: D c (ui) vs u> for the 
HCL with e c = 1.5, U ~ 7 (x = 3.5), r) = 0.8 and V 2 = 1.25, 
corresponding to n c ~ 0.2. The conduction band spectrum 
now shows a sharp resonance in the vicinity of the Fermi level 
uj = (cf Figs. 2, 3) as explained in text; po(u> — 0) = l/\/n is 
indicated by an arrow. 



In Fig. 5, comparison is made between an /-electron 
spectrum for the PAM for the HCL (solid line) and its 
counterpart for the single-impurity Anderson model (dash- 
ed line); the bare parameters chosen for illustration being 
1] = (symmetric /-level(s)), U ~ 1.3 (x = 0.5) and 
V 2 = 0.2, e c = 0.3 (corresponding to a hybridization 
strength, Eq.(4.21), of irA = (irV) 2 p (-e c ) ~ 1). The 
spectra, specifically ttAoD^(uj), are compared on the all 
scales level in the main figure (vs uj = w/i*); and on the 
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low-energy resonance scale in the inset. The main figure 
shows in addition the impurity spectrum for e c = (with 
the same parameters otherwise), which is of course the 
fully particle-hole symmetric AIM. The first point to note 
here is obvious: excepting the low-energy sector, the PAM 
and corresponding AIM spectra are qualitatively very sim- 
ilar on the 'all scales' level (which is not specific to the par- 
ticular parameter set employed). This general character- 
istic is in agreement with results from a numerical rcnor- 
malization group (NRG) study of the PAM and AIM [8], 
see e.g. Figs. 1,4, 6 of [8]. Note further, in comparison to the 
fully particle- hole symmetric AIM (dotted line), that for 
the moderate value of U chosen in Fig. 5 the asymmetry 
in the conduction band (e c ^ 0) shows up weakly in the 
Hubbard satellites, both in their intensities and maxima 
(which arc not quite symmetrically positioned about the 
Fermi level). 

It is naturally in the low-energy behaviour that the 
PAM and AIM spectra differ significantly, as evident in 
the inset to Fig. 5. In particular, for the PAM the lattice 
coherence generates a pseudogap above the Fermi level, 
as indeed expected from the quasiparticle behaviour of 
Df{tu), Eq.(3.11b) (albeit that the relative weakness of 
the pseudogap here reflects the moderate U considered, 
see §5.3). However even at low energies one does not learn 
much from comparison of PAM and AIM dynamics on an 
absolute scale, e.g. vs u) = u)/t*, and for a given set 
of bare parameters {e c , V 2 , U, r/}. For the parameters cho- 
sen in Fig. 5, it so happens that the low-energy scales for 
the PAM and AIM are very similar (with quasiparticle 
weights Z s=a 0.1). But that is not generically so: as dis- 
cussed in §5.2 the PAM lattice coherence scale and the 
Kondo scale for the AIM will in general be quite different 
for given bare parameters [8,23] (they are after all physi- 
cally distinct models), whence comparison of the two on an 
absolute scale is barely informative. What is required by 
contrast — particularly in strong coupling where there is a 
pristine separation between asymptotically vanishing low- 
energy scale(s) and non- universal scales such as Aq, t* = 1 
or U — is comparison of the universal scaling behaviour of 
the two models; in which the dependence of the respective 
low-energy scales on bare parameters is thereby eliminated 
and the underlying scaling behaviour exposed. This we 
believe is the most convincing (and possibly only) way to 
establish a connection between the high-energy scaling be- 
haviour of the PAM/KLM and underlying single-impurity 
physics. That key issue is considered in §5.3. 



5.2 Low-energy scale 

We first consider briefly how the low-energy coherence 
scale for the PAM, lu l = ZV 2 (Eq.(3.10)), is found within 
the LMA to depend on the bare/material parameters (e c , 
V 2 , U, rf) in the strong coupling Kondo lattice regime where 
rif — * 1; and how it compares to the Kondo scale for the 



corresponding AIM, lok 



7 \f' 



(with Zi mp denoting 



the quasiparticle weight for the single- impurity model). 
The scales lul and luk are indeed found to be exponen- 
tially small in strong coupling (as opposed to algebraically 




Fig. 5. All scales comparison of local /-electron spectra for the 
PAM and Anderson impurity model: -kAqD (w) vs oj = ojjt* 
(with A = -irV 2 Po (~e c )); for e c = 0.3, U ~ 1.3(ai = 0.5), V 2 = 
0.2 and n = 0. Solid line: PAM. Dashed line: AIM. Inset: 
comparison on the resonance scales. The impurity spectrum 
for e c — is also shown on the main figure (dotted line, the 
particle-hole symmetric AIM). 



small, such as arises using perturbation theory in U or 
variants thereof such as modified (iterated) perturbation 
theory [14-17]); leading thereby to the clean scale-separation 
that is a prerequisite to the scaling considerations of §5.3. 
The essential findings here, discussed below, agree with 
the NRG study of [8] and results obtained from the large- 
N/ slave boson mean- field (SBMF) approximation [23]: (a) 
That the scales u>l and uk are found to have the same 
exponential dependence on the bare parameters, but (b) 
depend very differently on the conduction band filling n c ; 
the lattice scale u>l being enhanced relative to its single- 
impurity counterpart asn c -> 1, but strongly diminished 
for low n c . 

The material dependence of the AIM Kondo scale wk{oc 
u> m ) arising within the LMA can be obtained analytically 
in strong coupling by direct analysis of the symmetry 
restoration condition Eq.(4.13). That was considered ex- 
plicitly in [27] for the case of a general impurity (with 

/-level asymmetry embodied as usual in r\ = 1 — 2-4p), 
but a symmetric host band. It is straightforward to ex- 
tend the analysis of [27] to include the ^-dependence of 
the hybridization A{u>) arising from an asymmetric host 
conduction band (which as anticipated in §4.5 docs not 
affect the final answer). Noting that the AIM hybridiza- 
tion strength is given by A n = TrV 2 po(—e c ) (Eq.(4.21)), 
this yields 



u>k oc cxp 



U 



(l- v 2 ) 



SV 2 p Q (-e c ) f( V 2 ) 



(5.1) 



(with the proportionality determined simply by a high- 
energy cutoff [27]); where f{rf) = \[l + ^(I-t? 2 )] (<G 
[1, |] for asymmetries \rj\ € [0,1] relevant to the Kondo 
regime). The exchange coupling J for the corresponding 
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Kondo model, obtained from the AIM in strong coupling 
by a Schricffer- Wolff transformation [2], is given by J — 
V 2 [jtj\ + Jj^j\\ = W 2 /{U(1 - r, 2 )]; whence Eq.(5.1) is 
equivalently lok oc exp(— l/[2po( — e c)Jf])- As pointed out 
in [27] the exponent here differs in general by the factor 
/(?7 2 ) from the exact result for the Kondo model, being 
as such exact only for the symmetric case where n = 
(although note that / is slowly varying in rj, lying e.g. 
within 10% of unity for n < 0.6). That notwithstanding 
we regard recovery of an exponentially small Kondo scale, 
close to the exact result in an obvious sense, as non-trivial; 
and add that provided the scale is indeed exponentially 
small, so that a clean separation of low (universal) and 
non-universal energies arises, its precise dependence on 
the bare parameters is in essence irrelevant to the issue 
of scaling in terms of lo/lok (as seen in [27] for the AIM 
itself). 



J 10" 




U/V 



Fig. 6. PAM coherence scale ujl on a log-scale vs U/V 2 ; 
for the BL with n — 0, e c = 0.1 (solid circles) and e c = 0.6 
(solid squares). Comparison is made to corresponding results 
for the AIM Kondo scale ujk'- e c = 0.1 (open circles) and e c — 
0.6 (open squares). u>l is enhanced relative to ujk for e c = 
0.1 but diminished for e c = 0.6 (see also Fig. 7). The straight 
lines indicate the exponents given by Eq.(5.1), found for both 
models. Inset: same results now shown vs U/Ao (with Aq — 
irV 2 po(— e c ))- The strong coupling gradients are then common 
for different e c , as implied by Eq.(5.1). 



For the PAM the coherence scale ujl is likewise ob- 
tained from solution of the symmetry restoration condi- 
tion Eq.(4.13), in this case numerically following the pro- 
cedure detailed in §4.4 (and with wj, found to be propor- 
tional to the spin-flip scale ui m as noted in §4.4). In the 
strong coupling Kondo lattice regime wi is again found to 
be exponentially small, with its exponential dependence 
on the bare parameters the same as ujk for the AIM. This 
is illustrated in Fig. 6 (for the BL) where, with r\ = and 
for e c = 0.1 and 0.6, the resultant lul is plotted on a log- 
arithmic scale vs U/V 2 ; and compared to the counterpart 
lok results for the AIM itself. For given e c the asymp- 



totic PAM and AIM curves are parallel, indeed indicating 
common [//^-dependence for the exponents of the two 
scales. When plotted vs U/V 2 as in the main figure, the 
slopes for different e c clearly differ; but when shown vs 
U/Aq = U/(irV 2 po(—e c )) as in the inset to Fig. 6, the gra- 
dients for different e c are now common in strong coupling, 
as implied by the exponential dependence of Eq.(5.1). The 
dependence of the exponents on the /-level asymmetry, as 
in Eq.(5.1), may likewise be verified by varying r\. And the 
same exponential dependence, Eq.(5.1), is found whether 
the BL or HCL is considered. 




Fig. 7. PAM coherence scale u>l vs n c , for the BL with 
U/V 2 = 23 and n — (solid circles), compared to the AIM 
Kondo scale ujk (open circles). Inset: ujl/ujk vs n c . 



While the exponents of the scales u>l and u>k have the 
same dependence on bare parameters, their dependence on 
the conduction band filling n c (which itself is determined 
solely by e c in strong coupling, see Eq.(3.13)) is quite dis- 
tinct for the two models. This is evident already in Fig. 6 
but seen more clearly in Fig. 7 where, for U/V 2 = 23 (and 
rj = 0) we show u>l and u>k vs n c (their ratio being shown 
in the inset). For n c — > 1, lul/luk = F(n c ) > 1 and the lat- 
tice scale is enhanced over its AIM counterpart; while with 
decreasing e c (and hence n c ) ujl diminishes progressively 
in comparison to lok, such that F{n c ) — > as n c — > 0. As 
noted above this general behaviour is in agreement with 
NRG [8] and SBMF results [23]. It is by contrast quite 
distinct from results arising from a Gutzwiller variational 
treatment [24] in which the lattice scale is always enhanced 
over its AIM counterpart; or from approaches based on lat- 
tice extensions of the non-crossing approximation (NCA) 
[18,19] in which the lattice scale, while in general moder- 
ately enhanced compared to u>k, is essentially equivalent 
to the AIM scale. 



5.3 Scaling 

The preceding discussion of how u>l and ujk for the two 
different models depend on bare parameters has been in- 
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Fig. 8. Universal spectral scaling. D c (u>) and nAoD-' (u>) vs uj' — uj/lul for the HCL with e c = 0.2 and r] = 0. The spectra for 
U ~ 5.1 (dashed line) and U ~ 6.6 (solid line) collapse to a common scaling form as a function of uj/uil- Inset: corresponding 
/-spectra shown by contrast on an absolute scale, vs ui/t*. The main figures also show the asymptotic low-|u/| quasiparticle 
spectra (Eqs.(3.11)), dotted lines; to which the full scaling spectra reduce for |u/| < 1. 



eluded in part because of the interest it has hitherto at- 
tracted in the literature. As noted earlier however we re- 
gard this matter as quite subsidiary in comparison to the 
strong coupling scaling behaviour of the lattice model it- 
self, as now considered. 

First we show that in the Kondo lattice regime the 
LMA indeed leads to universal scaling of single-particle 
dynamics in terms of lv/u>l, for fixed e c and r\. This is illus- 
trated for the hypercubic lattice in Fig. 8, for e c = 0.2 and 
rj = 0. The inset to the figure shows the /-electron spec- 
trum irAoD* (u>) on an absolute scale, i.e. vs ui (= w/t*), 
for V 2 = 0.2 and two different interaction strengths, U ~ 
5.1 (x = 2.5) and U ~ 6.6 (x = 3.25). In either case 
the resultant conduction band filling n c ~ 0.78, in agree- 
ment with the asymptotic result Eq.(3.13) which shows 
that n c is determined by e c alone; and likewise rif ~ 0.99. 
As seen from the inset the two spectra are quite distinct 
on an absolute scale and dependent on the bare model 
parameters, reflecting the exponential diminution of the 
coherence scale wi with increasing U/V 2 as in §5.2 above. 
The main part of Fig. 8 by contrast shows both ttAqD^uj) 
and the corresponding conduction spectrum D c (cu), now 
vs lu' = lu/lul, from which collapse to common scaling 
forms and hence universality is clear. While scaling has 
been demonstrated here by considering fixed V 2 upon in- 
creasing U in the KL regime, it is as expected dependent 
solely on the ratio U/V 2 (the same scaling spectrum aris- 
ing for fixed U upon decreasing V 2 ). 

As discussed in §3, adiabatic continuity to the non- 
interacting limit requires that for sufficiently low u/ the 
scaling spectra should reduce to the quasiparticle forms 
Eq.(3.11). That this behaviour is indeed recovered cor- 
rectly by the LMA is also seen in Fig. 8, where the re- 
sultant quasiparticle spectra are shown for comparison 



(dotted lines, as given by Eq.(3.11) with e 



/ 



tf 



(e c ) ob- 



tained from Eq.(3.12) with rif = 1): in the vicinity of the 
Fermi level (u/ = 0), and up to \ui'\ ~ 1 or so, agreement 
with the quasiparticle behaviour is essentially perfect. For 
larger |u/| by contrast, an evident departure from this sim- 
ple low-o/ asymptotic behaviour sets in; in particular, the 
quasiparticle /-electron spectra for large \u)'\ are seen to 
decay much more rapidly (~ l/|u/| 2 ) than the full LMA 
results, which show instead slowly decaying spectral tails. 
The latter, which as shown below decay logarithmically 
slowly, are a key feature of dynamics (see Figs. 11-13); 
reflecting genuine many-body scattering/lifetime effects, 
setting in for \u>'\ > 1 — 10 and dominating the scaling 
spectra (as well as transport properties on correspond- 
ing temperature scales, see e.g. [34,35]). Here we note 
in passing that scaling spectra arising from a SBMF ap- 
proximation are just the quasiparticle forms themselves, 
and are evidently deficient except for the lowest energy 
scales; and similarly that dynamics arising from modi- 
fied (iterated) perturbation theory [16,17] amount to lit- 
tle more than quasiparticle form, and similarly lack non- 
trivial high-energy scaling behaviour [35]. We also add 
that the spectral substructure seen in Fig. 8 just above the 
upper edge of the gap is not a numerical artefact, and us- 
ing the LMA self-energies can in fact be understood physi- 
cally in terms of correlated 'strings' of /-electron spin-flips 
on distinct lattice sites. It is however destroyed thermally 
on temperature scales which are a small fraction of wj, 
itself (as will be shown in subsequent work), and as such 
is but a minor feature of dynamics that we do not pursue 
further here. 

The spectra shown in Fig. 8 display an evident gap ly- 
ing slightly above the Fermi level (strictly a pseudogap for 
the HCL), as found also in approaches based on lattice ex- 
tensions of the NCA [18-20]. That such behaviour arises 
is to be expected, for it occurs likewise in the quasiparti- 
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cle spectra Eqs.(3.11) (see also the discussion at the end 
of §2.1). Note however that this gap become 'fully devel- 
oped' only in the strong coupling Kondo lattice regime; for 
weaker interaction strengths outside the scaling spectrum 
it is by contrast incompletely formed and evident only as a 
weaker pseudogap, as seen clearly e.g. in Fig. 5. But for suf- 
ficiently strong coupling we find that a well developed gap 
always arises (as the quasiparticle forms would suggest). 
Such behaviour is also found in recent NRG calculations 
[8] for n c ~ 1, but not for significantly lower conduction 
band fillings - see e.g. Fig. 7 of [8] for D c (ut) with n c = 0.6 
where, by contrast, only weaker pseudogap behaviour is 
evident. However the spectrum e.g. in Fig. 7 of [8] is clearly 
not close to strong coupling behaviour; as evidenced both 
from the fact that the D c (uj) shown there departs signif- 
icantly from the free conduction band envelope well into 
non- universal energy scales £)(£*), and because the quoted 
n c = 0.6 is far from its asymptotic strong coupling value 
of n c = 0.48 (from Eq.(3.13) above) for the bare parame- 
ters specified. Further resolution of this matter is clearly 
required, but we suspect that the parameter regime con- 
sidered e.g. in Figs. 6, 7 of [8] was not sufficiently strong 
coupling to uncover a well developed spectral gap. 

The scaling spectra shown in Fig. 8 refer specifically to 
'LMAF as detailed in §4.3, on which we focus in this paper. 
In Fig. 9 however we compare the resultant /-electron scal- 
ing spectra itAqD$ (uj) with those arising from 'LMAIF, 
where (see §4.3) the polarization propagators entering the 
LMA self-energies are further renormalized in terms of 
the host/medium propagators {Q a (u>)}; again for the HCL 
with e c = 0.2 and r] = 0. The inset shows the LMAI/II 
comparison out to a/ = 10, while the main figure extends 
to much larger scales. And the two spectra are seen to be 
essentially coincident on all u)' scales (as they ought to be 
if the LMA captures adequately the scaling spectrum). 




-25 

oo/ox 



Fig. 9. Comparison of /-electron scaling spectra arising from 
LMAI (solid line) and LMAII (dashed line), as explained in 
text: ttAqD (w) vs cj/lul for e c = 0.2 and n — 0. The two 
levels of LMA yield essentially coincident scaling spectra. 



Fig. 8 above illustrates that universal spectral scaling, 
independent of U and V 2 , arises for fixed e c and r\ which 
embody respectively asymmetry in the conduction band 
and /-levels. This we find to be quite general: D c (uj) and/or 
7rZ\o-D^(w) exhibit scaling as an entire function of u/ = 
lu/u>l only for fixed (e c ,rj), i.e. distinct scaling spectra 
arise for different {e Cl rj). Much more subtly however, the 
e c - and //-dependences of the scaling spectra depend upon 
the w'-regimes considered, as now explained. We begin 
with the simple case of low-a/. As pointed out in §3, the 
quasiparticle spectra Eq.(3.11) imply that — in their w'- 
regime of validity — the scaling spectra should actually 
be independent of the /-level asymmetry r\. That this is 
recovered within the LMA is seen in Fig. 10 for the HCL 
where, for fixed e c = 0.3, the /-electron scaling spectra 
ttAqD^ (uj) arc shown for three different /-asymmetries 
i] = 0, 0.2, 0.4. For \u>'\ < 1, the regime where the quasipar- 
ticle forms hold, the LMA scaling spectra are indeed seen 
to be independent of rj; while for larger- |u/| by contrast, 
an ^-dependence to the spectra is evident (and discussed 
further below). Likewise, for \u'\ < 1, the quasiparticle 
forms Eq.(3.11) show that the scaling spectra depend ex- 
plicitly on e c , as well as on the underlying lattice itself 
(embodied in the specific form for po(u>)). 




Fig. 10. Independence of the low-|a/| scaling spectra on 
the /-level asymmetry, n. For the HCL with fixed e c — 0.3, 
ttAoD-' (lu) vs w' = co/i-ol is shown for rj — (solid line), 0.2 
(dashed line) and 0.4 (point dash). The limiting low-|u>'| quasi- 
particle form is again shown for comparison (dotted line). 



But what of higher energies in the scaling spectra? 
Here as we now show the low-u/ situation above is re- 
versed: the high-energy scaling behaviour of the /-electron 
spectrum is dependent on the asymmetry 77, but indepen- 
dent of both e c and the underlying host lattice; the latter 
in turn being intimately related to the emergence of ef- 
fective single-impurity physics in the high-energy scaling 
behaviour of the PAM. 

To see this, Fig. 11 (for the HCL) shows wA D f (u) vs 
lu' = lu/lul up to \u>'\ = 500, for 77 = and three dif- 



18 



N.S.Vidhyadhiraja and David. E.Logan: Dynamics and scaling in the periodic Anderson model 



0.8 



3, 0.6- 

I— I 

<°o.4 



0.2 





I i I i I 

-100 -50 50 100 



-400 



-200 




00/GX 



200 



400 



Fig. 11. TvA D f (w) vs J = lu/ljl for the HCL up to \u'\ = 
500, for r) = and three different e c — (solid line), 0.3 (dashed 
line) and 0.5 (point dash); obtained explicitly for U ~ 6.6 and 
V 2 = 0.2. Inset: irA D f (u>) vs J for fixed e c = 0.3 with 
increasing interaction strength x = ^£/|/u| — 2.0 (solid line), 
2.5 (double point dash) and 3.25 (long dash); the true scaling 
limit is also indicated (dotted line). Full discussion in text. 



ferent e c — 0,0.3,0.5 (progressively diminishing conduc- 
tion band filling n c ); the particular results shown having 
been obtained explicitly for U ~ 6.6 (a: = \U\[i\ = 3.25) 
and V 2 = 0.2. Looking at the negative-u/ side in partic- 
ular, it is clear that the slowly decaying spectral 'tails' 
are indeed asymptotically common for the different e c 's. 
On the positive-u/ side there might appear from the fig- 
ure to be a residual weak dependence of the spectral tails 
on e c . That however is simply a reflection of the natu- 
ral fact that the value of U/V 2 required to reach the full 
asymptotic scaling spectra is dependent upon e c . This is 
illustrated further in the inset to Fig. 11, which for fixed 
e c = 0.3 shows (on an expanded scale) the evolution of 
the scaling spectrum with increasing interaction strength: 
x = 2.0,2.5 and 3.25. Looking at the positive u/ side one 
sees that the true scaling limit (dotted line) is steadily ap- 
proached upon increasing the interaction strength, but not 
reached until a U somewhat in excess of ~ 6.6 (x = 3.25). 
On the negative lu' side by contrast, the scaling limit is 
already reached by a; = 2.5 (U ~ 5.1) and does not change 
with further increasing interaction strength. This is why 
the e c -indcpcndcncc of the spectral tails is clearly evident 
only on the u/ < side of the main figure; upon increas- 
ing U however both sides of the scaling spectra show this 
behaviour. 

What then is the functional form of the large- |u/| spec- 
tral tails? On physical grounds one expects that on suf- 
ficiently high energy and/or temperature scales, the /- 
electrons in the Kondo lattice regime of the PAM should 
be screened in an essentially incoherent single-impurity 
fashion; and thus that effective single-impurity physics 
should arise in the lattice model at high energies, quite dis- 
tinct from the effects of lattice coherence evident on low- 



energy scales lo/lol ~ 0(1). For the AIM itself the spectral 
tails of the local impurity scaling spectrum Di mp {oj) can 
be obtained analytically within the LMA [27,28]; being 
given by 



irA D imp (u>) ~ - 



JMc|w| 



1 LlMc 



f 25 
(5.2) 



(shown explicitly for r\ = 0) where u> = lu/lok with luk the 
Kondo scale for the AIM discussed in §5.2, and c a pure 
constant 0(1). For w > 5-10 or so, Eq.(5.2) is known [28] 
to describe quantitatively the spectral tails arising from 
NRG calculations; the exact high energy scaling asymp- 
tote irAoDi m p(uj) ~ 37r 2 /[161n (lu/ojk)} being recovered 
in particular. If effective single-impurity behaviour arises 
in the PAM on energy scales encompassed by the en- 
sealing regime, then the spectral tails thereof should have 
the same scaling form as for the AIM, i.e. should be given 

by 



wA D f (u 



1 



;|ln(a|^|)]2 + l [lln(o|w'|)] 2 + 25 



(5.3) 



with u>' = oj/u>l and a a constant 0(1). Note that such 
comparison requires neither a knowledge of how the low- 
energy scales for the two distinct models (wl and uik) de- 
pend on the bare material parameters, nor any assumption 
that the Kondo scale for the AIM itself is at all relevant 
to the PAM; points to which we return again in §5.4. 
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Fig. 12. Effective single-impurity physics arising in the PAM 
scaling spectra at high energies, as explained in the text. The 
PAM -K AqD* (uj) vs \lo\/lol on an expanded vertical scale, for 
n — and e c — (solid line), 0.3 (dashed line) and 0.5 (double 
point dash); compared to the scaling form Eq.(5.3) (dotted 
line, barely distinguishable in the figure). 



Eq.(5.3) indeed describes the tail behaviour of the PAM 
scaling spectra, and as such establishes the connection to 
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effective single-impurity behaviour at high energies. This 
is shown explicitly in Fig. 12 where, again for n = and 
e c = 0,0.3 and 0.5 (as in Fig. 11), ttAqD '(co) is shown 
vs lu' on an expanded vertical scale, and compared to 
Eq.(5.3) (with the constant a ~ 0.55 determined numer- 
ically). That form is clearly seen to hold asymptotically 
for the different e c 's; and even for lower lu' the spectra 
display only a very weak dependence on e c . Neither is it 
lattice dependent, the same asymptotic tail behaviour be- 
ing found to arise whether the HCL or BL is considered; 
naturally so, for effective single-impurity scaling physics 
should be independent of the 'host' lattice. 




Fig. 13. /-electron scaling spectra itAqD (u>) vs lJ = w/wl, 
for n = (solid line) and r\ = 0.3 (dashed line); shown explicitly 
for tc — (the scaling spectra at low- a/ depend upon e c , see 
text). 



While the discussion above has focused on varying e c 
(and hence n c ) for symmetric /-levels r] = 0, the behaviour 
found is quite general. For fixed n ^ the high-energy 
PAM scaling spectrum is likewise independent of both e c 
and the lattice type; and is again found to have precisely 
the same scaling form as its counterpart for the AIM (the 
generalisation of Eq.(5.2) to finite-??, specifically Eq.(5.5) 
of [27]). As for the AIM [27] the resultant /-electron scal- 
ing spectra arc now //-dependent, as illustrated in Fig. 13 
where ttA D^ '(lu) is compared for r\ = and 0.3 (and is 
shown specifically for e c = 0, bearing in mind that the 
spectra at low-o/ depend on e c as discussed above). The 
77-dependence of the spectral tails is clearly evident, albeit 
rather weakly so for positive u)' in particular. 

The above results capture the evolution of the scal- 
ing spectra appropriate to the Kondo lattice regime of 
the PAM, from the low-energy behaviour symptomatic of 
the coherent Fermi liquid state through to the effective 
incoherent single-impurity physics found to arise at high 
energies — but still in the lu' — lu/lul scaling regime. Fi- 
nally, we add that while our exclusive focus here has been 
on single-particle dynamics, the results obtained naturally 
have direct implications for transport and optical proper- 



ties of heavy fermions; these will be considered in a sub- 
sequent paper. 

5.4 Discussion: Exhaustion? 

The discussion of the previous section brings us to Nozieres' 
issue of 'exhaustion' [37,38] and the question of how a co- 
herent Fermi liquid state forms in a concentrated Ander- 
son/Kondo lattice. For a single impurity Anderson/Kondo 
model, assuming [38] the only conduction electrons eli- 
gible to provide Kondo screening are those lying within 
~ u>k (= '2~W') of the Fermi level (lu = 0), the number 
of such is Ns ~ NlcIq(0)u)k', with Nl the total num- 
ber of lattice sites (and (Iq(cu) = po(uj — e c ) the free con- 
duction band dos, normalised to unity). In the strong 
coupling Kondo regime, cIq(0)u)k is of course exponen- 
tially small; but Ns, the number of available screening 
electrons per the single impurity spin, obviously remains 
macroscopically large. That situation changes drastically 
in the concentrated Anderson/Kondo lattice. Now there 
are Nl spins (/-electrons) to screen; so the number of 
electrons per f-spin available to provide Kondo screen- 
ing is Ns/Nl <~ cIq(0)ujk — itself exponentially small. 
That raises the issue of 'exhaustion' [37,38]: how so few 
screening electrons lead to the formation of a coherent 
Fermi liquid ground state. Nozieres' argument [38] is that 
this effectively arises through a two-stage process with de- 
creasing energy/temperature scale. Neglecting the RKKY 
interaction, on high energy/temperature scales the local 
/-spins are first Kondo screened in an essentially incoher- 
ent, single- impurity fashion; while with further decreasing 
energy this effective single-impurity regime crosses over 
into lattice coherent behaviour through collective screen- 
ing/isotropization of the / spins. Two underlying scales 
are then argued to emerge: a high energy single-impurity 
Kondo scale luk corresponding to the incoherent effec- 
tive single-impurity physics; and a second, lower lattice 
scale lol (= 'T c ') signifying the onset of lattice coherence. 
Nozieres has provided intuitive arguments [38] to suggest 
that, at most, lul ~ d%(Q)[tu k] 2 ', which, since luk itself is 
exponentially small, means that lul and luk are radically 
distinct in scaling terms (as elaborated below). Further, 
as noted by Pruschke et. al. [8], Nozieres' phenomcnologi- 
cal arguments are not in fact particular to low conduction 
band filling n c , and if correct imply two-scale exhaustion 
physics should be the generic situation. 

Much work has since ensued on the question of exhaus- 
tion [8,11,12,17,23] via DMFT studies of the PAM and/or 
KLM. Regarding scales per se there appears now to be an 
approaching consensus [8,23] that lul oc [luk] 2 does not 
arise; but rather that u>l oc u>k, with a proportionality de- 
pendent on the conduction band filling n c (or equivalcntly 
e c ): lul/luk = F(n c ), with which the present work concurs 
as in §5.2. That granted however, it has nonetheless still 
been suggested e.g. in [8,23] that away from half-filling 
where F(n c ) < 1, a two-scale picture arises. This under- 
lies the qualitative notion of 'protracted screening' from 
Quantum Monte Carlo/Maximum Entropy Method stud- 
ies [11,12], where the existence of two scales is inferred 
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because the thermal evolution of dynamics is slower (or 
'protracted') for the PAM than the AIM; and it arises 
likewise in the large-Af mean-field study of [23] (where 
the emphasis is on low conduction band filling). It is also 
suggested in the NRG study of [8] (§V of [8]), although we 
add that the NRG results themselves were not argued in 
[8] to provide evidence for relevance of the single-impurity 
scale lok to dynamics of the PAM itself. 

As is evident from the results of the previous section, 
we dissent from the view of a two-scale picture (without 
disagreeing with specific results obtained e.g. in [8,23]). 
That deserves a careful explanation. The first point to 
make here is that if if two distinct scales were relevant to 
the PAM/KLM, in particular to its behaviour all the way 
from the coherent Fermi liquid through to effective single- 
impurity physics, then there should be two distinct scaling 
regimes of the model. But that could only be ascertained 
by investigation of the scaling behaviour of physical prop- 
erties (in the relevant strongly correlated regime); which 
has not hitherto been considered. 

Let us then suppose that two distinct scales arise. The 
lower will be the coherence scale u>l, and (without preju- 
dice as to its origin) call the higher energy scale loh ■ As a 
function of the bare model parameters it is taken as read 
that both scales vanish asymptotically in strong coupling 
{i.e. become exponentially small, well separated from non- 
universal scales). Now suppose the ratio ujh/^l — * oo 
in strong coupling (which would be the case e.g. in the 
Nozieres exhaustion scenario [38] where luh corresponds 
to the AIM scale lok and u>l oc [oj/<-] 2 ). In that case, scal- 
ing of physical properties in terms of u/ = oj/ujl (or T/ojl) 
would project out to infinity energies on the scale of u>h 
(as well as the usual irrelevant non- universal scales); this 
is the ''lu/ojl scaling regime. By contrast, scaling in terms 
of to" = lu/ojh would project out both non-universal scales 
(to infinity as usual) and all energies on the scale of ujl or 
finite multiples thereof (to zero); this is the second, 'w/wh' 
scaling regime. The coherent Fermi liquid will of course be 
encapsulated in the former, lu/lul scaling regime. Regard- 
ing the crossover to effective single-impurity behaviour in 
the PAM/KLM, two obvious possibilities arise, (a) That 
this crossover is set by the ojh scale (as in the Nozieres pic- 
ture). It will then arise only in the second scaling regime; 
effective single impurity behaviour will not thus be evi- 
dent as a function of uj/lul- (b) That the crossover occurs 
in the lo/lol scaling regime. In that case the second puta- 
tive scaling regime is irrelevant — at least to the central 
issue of understanding the evolution of the PAM/KLM 
from the coherent Fermi liquid through to effective single 
impurity behaviour. 

The behaviour hypothesised above reflects a genuine 
two-scale description. But if by contrast the ratio luh/i^l 
is a constant in strong coupling then the scales u>l and 
<jJh are equivalent, differing only by the constant but fun- 
damentally equivalent in scaling terms. In this case there 
is no essential distinction between ojl and loh, and obvi- 
ously only one scaling regime. Such behaviour is of course 
well known to arise e.g. in the Anderson impurity model, 
where the 'Kondo scale' appears in many different but 



equivalent guises: e.g. [2] the usual Kondo temperature 
Tk obtained from the impurity susceptibility, the half- 
width at half maximum of the Kondo/Abrikosov-Suhl res- 
onance in the single-particle spectrum, or AoZi mp with 
Zimp the impurity quasiparticle weight. Each is propor- 
tional to the other; all are manifestations of the single 
underlying Kondo scale. 

Within the present theory the results of §5.3 show 
that a one-scale picture arises: the evolution from a coher- 
ent Fermi liquid to effective single-impurity physics arises 
clearly when the scaling behaviour of dynamics is con- 
sidered as a function of uj/lol (and we find no evidence 
for a 'higher' scaling regime). Neither is such behaviour 
confined to single-particle dynamics. As will be shown 
in a subsequent paper the resistivity p(T), including its 
crossover to effective single-impurity behaviour, likewise 
exhibits one-parameter universal scaling in terms of T/lul. 
In parallel to the above comments on the AIM, that does 
not of course preclude the coherence scale appearing in 
different but equivalent guises. For example the peak max- 
imum in the resistivity, often chosen as a measure of the 
low-energy scale, is not identically u)l\ but the two are 
simply proportional, indicative of one-parameter scaling. 

One further point should be noted here. The connec- 
tion to effective single-impurity physics in the PAM/KLM 
established in §5.3 arises from comparison of the scaling 
forms of the spectra for the PAM/KLM and the AIM, i. e. 
as functions of lj/u>l and uj/lok respectively. That does 
not require any knowledge of how the separate scales for 
the two models, u>l and u>k, themselves depend on the 
bare material parameters. Neither docs it require any as- 
sumption that the Kondo scale luk for the single-impurity 
model itself is directly relevant to the PAM/KLM (the 
two are after all different models). It is nonetheless the 
case that, from specific study (§5.2) of how the respective 
scales for the two models depend upon the bare parame- 
ters, we find ojl/luk = F(n c ) in agreement with previous 
work [8,23]; and while that behaviour is neither required 
nor relevant in establishing via scaling the connection to 
effective single-impurity physics in the PAM/KLM, it does 
mean that the two scales are fundamentally equivalent if 
used formally as scaling variables, i.e. one may choose 
equivalently to use lok to scale the PAM spectra. Hence, 
while we concur with the basic results of [8,23] for the re- 
lation between ujk and u>l , we naturally disagree with the 
view that a two-scale picture arises. 

We also note parenthetically that the single-scale pic- 
ture found here simply obviates an apparent conundrum 
raised in [8] , namely how one rationalises the regime n c w 
I where lul/ujk > 1 — the lattice coherence scale is now 
larger than the single-impurity Kondo scale. As pointed 
out in [8] this presents a self-evident problem for interpre- 
tations based on the assumption that the AIM u>k is the 
relevant scale for effective single-impurity physics in the 
PAM/KLM, while wi sets the scale for lattice coherence. 
From the viewpoint of the present work however, this is 
obviously not an issue; effective single-impurity physics in 
the PAM/KLM arises as naturally for n c — > f (see also 
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Fig. 14. Weak coupling LMA /-electron spectrum (solid line) vs co/t, for the Bethe lattice, with U — 0.75 (right panel) and 
U — 0.25 (left panel); and for e c — 0.2, r\ = and V 2 = 0. Corresponding results from second order perturbation theory in U 
are also shown (dashed lines); for the lower-f/ shown this is indistinguishable from the LMA spectrum. 



[34,35]) as it docs for lower conduction band fillings n c 
where u>l/ujk < 1. 

The obvious conclusion from the preceding discussion 
is that we find no compelling evidence for two-scale ex- 
haustion. A sceptic can naturally argue that since the 
present theory is approximate it is open to doubt. That is 
of course true, as it is for any theory. But the evidence here 
certainly points away from exhaustion and, should further 
support for the idea be forthcoming, it will in our view re- 
quire convincing scaling arguments to be established. 



6 Concluding remarks 

We have developed in this paper a local moment approach 
to single-particle dynamics of the periodic Anderson model 
within the framework of dynamical mean-field theory, for 
the generic asymmetric case relevant to heavy fermion 
metals. For obvious physical reasons our primary inter- 
est has been the strongly correlated Kondo lattice regime, 
of essentially localized spins n/ — ► 1 but with general 
conduction band filling n c , which the intrinsically non- 
perturbative nature of the LMA renders readily accessi- 
ble. The exponentially small lattice coherence scale lot, 
inherent to the Kondo lattice regime leads in particular 
to a clean separation of low-energy ('universal') and high- 
energy scales, and hence to universal scaling behaviour 
of dynamics. This has been a central focus of the present 
work and a rich description of scaling spectra results, span- 
ning all lu /lo L-scsdes. With increasing u/ = lj/lul, dynam- 
ics are found to cross over from the low-energy quasipar- 
ticle behaviour symptomatic of the coherent Fermi liq- 
uid state to essentially incoherent single-impurity Andcr- 
son/Kondo scaling physics at high-a/. The former, low-u/ 
behaviour depends naturally on both the conduction band 
filling and underlying lattice 'type'. The latter by con- 
trast depends on neither, consistent with one's physical 
expectation of effective single-impurity physics; and the 



crossover from coherent Fermi liquid to effective single- 
impurity behaviour in the PAM, established as it is by 
scaling, neither presumes nor requires any particular rela- 
tion between the PAM coherence scale, u>l, and the Kondo 
scale lok for the corresponding 'real' AIM arising when 
only a single /-level is coupled to the conduction band. 

While our almost exclusive emphasis has been on the 
strongly correlated regime we add that, as for the Ander- 
son impurity models considered hitherto [26,27], all inter- 
action strengths are nonetheless encompassed by the LMA 
including simple pcrturbativc behaviour in weak coupling 
- an illustration of the latter being given in Fig. 14 where, 
for U = 0.75 and 0.25 (with e c = 0.2, r\ = and V 2 = 0.2) 
LMA results for the /-electron spectrum ttAqD^lu) vs 
u>/t* are compared to those arising from second order per- 
turbation theory (PT) in the interaction U: with decreas- 
ing interaction the LMA spectrum clearly reduces to that 
arising from PT, being indistinguishable from it for the 
lower U shown. 

The essential criteria for a successful description of the 
PAM thus appear to be met by the LMA; all energy scales, 
and interaction strengths from weak to strong coupling, 
being handled in a unified theoretical framework. Further 
(as mentioned in §1), a description of the single-particle 
dynamics considered here is a prerequisite to determining 
transport and optical properties of heavy fermions; which 
subject, including explicit comparison to experiment, will 
be considered in a forthcoming paper. 



We express our thanks to the EPSRC for supporting this re- 
search. 
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